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ABSTRACT 

Using the Infrared Spectrograph on the Spitzer Space Telescope, we observed the Antennae galaxies 
obtaining spectral maps of the entire central region and high signal-to-noise 5 — 38/im spectra of the 
two galactic nuclei and six infrared-luminous regions. 

The total infrared luminosity of our six IR peaks plus the two nuclei is Ljr — 3.8 • 10^*^ L©, with their 
derived star formation rates ranging between 0.2 and 2 Mo/yr, with a total of 6.6M0/yr. None of the 
typical mid-infrared tracers of AGN activity is detected in either nucleus of the system, excluding the 
presence of an dust enshrouded accretion disk. The hardest and most luminous radiation originates 
from two compact clusters in the southern part of the overlap region, which also have the highest dust 
temperatures. PAH emission and other tracers of softer radiation are spatially extended throughout 
and beyond the overlap region, but regions with harder and intenser radiation field show a reduced 
PAH strength. 

The strong H2 emission is rather confined around the nucleus of NGC 4039, where shocks appear to 
be the dominant excitation mechanism, and the southern part of the overlap region, where it traces 
the most recent starburst activity. The luminosity ratio between the warm molecular gas (traced by 
the H2 lines) and the total far-IR emission is ~ 1.6 • 10""*, similar to that found in many starburst 
and ULIRGs. The total mass of warm H2 in the Antennae is 2.5 • 10^ M0, with a fraction of warm 
to total H2 gas mass of about 0.35%. The average warm H2 temperature is 302 ± 26K and appears 
anti-correlated with the radiation field hardness, possibly due to an evolution of the PDR morphology. 
The previously reported tight correlation between the H2 and PAH emission was not found but higher 
total PAH emission to continuum ratios were found in PDRs with warmer gas. 

Subject headings: Telescopes: Spitzer, Galaxies: interactions, starburst, star clusters; ISM: dust, 
extinction, HII regions, molecular hydrogen; Infrared: galaxies 



\ 1. INTRODUCTION 

The Antennae galaxies, NGC 4038/39 (Arp 244), are 
arguably the best-known example of interacting galaxies. 
The two galactic nuclei of NGC 4038 and NGC 4039, and 
the dust-obscured overlap region in between, exhibit one 
of the most stunning examples of starburst activity^ in 
the nearby Un iverse. They are the first system in the 
■ iToomrd (fl97l merger sequence. In the Antennae we 
observe the symmetric encounter between two normal Sc 
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^ Throughout the paper we use both t erms ULIRGs and star- 
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term 'starburst' for the high star formation activity in NGC 7714). 
ULIRGs, unless dominated by AGN activity, constitute the sub- 
group of the most lumin ous starbursts with a tota l infrared lumi- 
nosity of LiR > 10^2 L0 llSanders fc MirabelllTog^ ') . 



type spirals at an early merger stage (|Mihos &: Hernguistj 

[TgM). 

From the four IRAS bands iSanders et all (|2003D 
derived total infrared luminosities of logipiR = 
^40-400^™ — 10.73 and logLjR = is-ioooAim = 10-84, 
corresponding to Lfir = 5.6 • 10^° L0 and Lir — 7.2 • 
10^° Lq for a distance of 22 Mpc. With a luminosity 
just below 10^^ L0 the Antennae does not classify as a 
luminous infrared galaxy (LIRG). 

The distance to the Antennae has been a sub- 
ject of re cent cont roversy. Earlier estimates by 
I'Wilson et al. ( 2000): i Gao et al.l (|2001[ ) found 20 Mpc, 
and lWhitniore et al . (1999) determined 19 .2 Mpc. How- 
ever, using HST-ACS . Saviane et all ([200l observed the 
tip of the red giant branch of the stellar populations in 
the southern tail of the Antennae, and determined a dis- 
tance of only 13.3 ± 1.0 Mp c. This estimate was soon 
again altered when, recently, iSchweizer et al.l (120081 ) de- 
rived 22.3 ± 2.8 Mpc from the type la supernova 2007sr 
light curve, much closer to the ol d estimates, which were 
mainly based on the Hubble flow. ISchweizer et al.l (j2008D 
also re-analyzed the HST-ACS data and identified a dif- 
ferent location of the tip of the red giant branch, which 
agreed with their SNla distance. Combi ning three in- 
dependent methods ISchweizer et al.l (j2008f ) determined a 
distance of 22±3.0 Mpc, which is the distance we assume 
throughout this paper. 

The first deep optical images of the Antennae 
taken with the Wide Field Ca mera on HST revealed 
over 700 point-like objects (|Whitmore fc Schweizeil 
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I1995D . Subsequent V-band observations with WFPC2 
increased the sensitivity by three magnitudes and 
revealed between 800 and 80 00 clusters in four age 
ranges (|Whitmore et al.l I1999D . These and related 
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119991: iWhitmore & Zhand 120021: 


Kassin et al. 120031: 


Brandl et al.l 120051; Mengel et all 


20051 iFall et all 120051; iBastian et all I2006D started 



extensive studies of the properties of extragalac- 
tic star clusters in a statistically significant way. 
The Antennae provide an excellent environment to 
study the transformation of supergiant molecular 
clouds (SGMCs) into super star clusters (SSCs) 
()de La Fuente Marcos &: de La Fuente Marcoi l2006l) 



and their further dy namica l evolutio n (Whitmor e et al 
19991; IFall et all [2005; Bastian et all ,2006; Anders et al 



2007t Vhitmo reet all |2007|). 



ies (TCilbert et al 



Numerou s othe r stud- 

^_ iMeneel et all l200ll l200l : 

iGilbert fc Graham! 120071 ) focused on detailed investi- 
gations of a few individual SSCs via means of near-IR 
spectroscopy. 

The first high signal-to-noise mid-IR images of the 
Antennae were made by the Infrared Space Observa- 
tory ISO. ISOCAM observations at angular resolutions 
of 5 — 8" showed that the overlap region contributes 
more than half of th e total luminosity ob served in the 
12.5 — 18^m range (|Vigrouxet al.ll 19961 ). The com- 
parison between HST and ISOCAM images revealed a 
compact, optically-obscured knot which produces 15% o f 
the total 12.5 - 18^m luminosity (jMirabe l et al."1998l. 
Further ISO studies were p e rformed by Ku nz e et al. 
(Il996l); iFischer et al.l (119961 ): iKlaas et alf (|l997f ) and 
iHaas et al.l (I2005D followed by Spitzer-IRAC 3 - S/xm 
observations (jWang et al.ll2004 l. Recent sub-arcsecond 
imaging and N-band spectroscopy with VLT-VISIR re- 
solved the Hii region / PDR interface and revealed a 
highly obscured SSC that does not have an optical or 
near-IR counterpart (jSniiders et al.ll2006l ). 

Observations at the longer far-infrared, sub-millimeter 
and radio wavelengths generally agree with the mid-IR 
jicture. The early works of Hummel fc van der Hulst) 



(|1986[ ) with the VLA at 1.465 GHz and 4.885 GHz found 
thermal radio knots, which account for about 35% of 
the total emission, coinciding with peaks in H^ emis- 
sion associated with recent star form ation. The high 
resolu tion radio maps with the VLA (|Nefr fc UlvestadI 
l2000f) at 6 and 4 cm became t he ba sis for many sub- 
sequent studies. iNikola et al.l (jl998D mapped the An- 
tennae in the [C II] cooling line associated with photo- 
dissociation regions at an angular resolution of 55" and 
found that the starburst activity is confi ned to small 
regions of high star formation efficiency. 'Wils on et al.l 
(|2000l [2003.) found an excellent correlation between the 
strengths of the CO emission and the 15/im broad 
band emission seen by ISO, and determined masses of 
(3 — 6) • 10® M0 for the largest molecular complexes, typ- 
ically an order of magnitude larger than the largest struc- 
tures found in the disks of more quiescent spiral galax- 
ies. The total molecular gas mass^° in the overlap re- 

These estimates generally use the standard Galactic conver- 
sion factor of Mjj^ = 4.78 X (^Lqo / Kkms~^pc'^) Mq to derive 
the total H2 mass from the CO luminosity, multiplied by 1.36 to 
account for the heavier elements. 



gion is 1.2 ■ 10^ Mq (jStanford et al.lll990D - several times 
higher than the amount of gas in the two nuclei - while 
the molecular ga s mass for the e ntire Antennae system 
is 1.5 • IO^^Mq ()Gao et al.ll200lD. At a star formation 



rate (SFR) of ~ 20M,r. vr- ~(IZhang. Fall fc Whitmord 
20011). the gas c onsumption timescale is only 700 Myr. 



Gap et aTl (j2001f ) also estimated a "normal" star forma- 



tion efficiency (SFE) over the entire Antennae system 
of Lin/M/fj ^ 4.2L0/M0, which is simil a r to t hat of 
GMCs in the Galactic disk. ISchulz et all (|2007f) com- 
bined i^CO (1-0), (2-1), (3-2) and 870Aim maps with data 
from X-ray to radio and PDR models, and found that the 
clouds have dense (4 x 10'*cto~^) cores and low kinetic 
temperatures (< 25 K), showing no signs of intense star- 
burst activity. 

Finally, numerous X-ray observations, predominantly 
with Chandra, have revealed many X-ray sources, in- 
cluding several lon g-term variable, ultra-lu minous X- 
ray (ULX) sources (jZezas et al.ll2002ai [2006I ). Although 
most ULXs are likely black hole/high mass X-ray binaries 
accreting via Roche lobe overfiow, Feng fc Kaaret (200^ 
have argued that the most luminous ULX, X-16, is a can- 
didate for an interm ediate mass black hole (IMBH), while 
IZezas et"an ()2007D note that its X-ray luminosity could 
be produced by a ^ 80 Mq black hole accreting at the 
Eddington limit. From the spati al distribution of th e in- 
frared counterparts to the ULXs. lClark et all ()2007( ) con- 
cluded t hat they ma i nly trace the recent star formation 
history. iBaldi et al.l (|2006a) studied metal abundances 
in the hot interstellar medium and found variations from 
0.2 Z0 to 20 — 30 Zq, but found no correlation between 
the radio/optical star f ormation indicators and metallic- 
ity (|Baldi et al.ll2006bf ). 

In this paper we study the central region of the An- 
tennae galaxies based on spatially resolved mid-IR spec- 
troscopy. Our aim is to characterize the starburst activ- 
ity in general, the super star clusters in particular, and 
their interplay (heating and ionization) with the inter- 
stellar medium. From the mid-infrared spectra we derive 
a total infrared luminosity for each of our six star forming 
peaks. These luminosities are used to estimate the star 
formation rates. Particular attention will be given to the 
H2 rotational lines, which probe the warm component of 
molecular hydrogen and which will be used to derive the 
gas masses and temperatures. The outline of this paper 
is as follows: In the next section we describe the obser- 
vations along with the data reduction and calibration. 
Section |3| discusses how the spectral features have been 
measured and the key parameters were derived, and sec- 
tion m presents the results, both on individual clusters 
and on the central region as a whole, followed by the 
summary. 

2. OBSERVATIONS AND DATA REDUCTION 

We used the Infrared Spectrograph^^ {IRS) 
(iHou ck et al.ll2004D o n board the Spitzer Space Telescope 
(Wer ner et al.l 12004 ) to observe the central interaction 
region of the Antennae galaxies at unprecedented depth. 
The IRS observations of the Antennae are part of a 
guaranteed time program (PI Houck, Spitzer PID 21) on 

The IRS was a collaborative venture between Cornell Univer- 
sity and Ball Aerospace Corporation funded by NASA through the 
Jet Propulsion Laboratory and the Ames Research Center. 
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TABLE 1 
Observing Parameters 



'hires' 



'lores' 



Obs.date 



AOR'^-keys 



Slit sizes 



3 Jan 2005 



3842304 
3842560 



24 X 30s (SH) 
10 X 60s (LH) 
4."7 X ll."3 (SH) 
ll."l X 22."3 (LH) 



3 Jul 2005 
12 Jul 2005 
23 Jan 2006 
3841536 
3841792 
16707840 
16708096 
16708352 
2 X 14s (SL) 
2 X 14s (LL) 
3."6 X 57" (SL) 
10."5 X 168" (LL) 



Astronomical Observing Request*^ Number of expo- 
sures per position X exposure time 

the spectroscopic study of star formation in interacting 
galaxies. The observing parameters are Hsted in Table [1] 
and described in more detail in Sections 12.11 and 12.21 




2.1. 'Hires' Observations 

We have made observations in IRS high resolution 
('hires') mode at i? = 600 with the Short-High (SH) 
[9.9 - 19.6Mm] and Long-High (LH) [18.7 - 37.2/im] 
modules. The angular slit sizes are listed Table [1] and 
correspond to approximately 500 x 1200 pc^ (SH) and 
1200 X 2380 pc^ (LH) at the distance of the Antennae. 

The spectra from the IRS 'hires' modules were taken 
with a standard staring mode Astronomical Observing 
Template (AOT), producing two exposures per 'cycle' at 
separate nod positions along the slit. The coordinates 
of the eight observed regions within NGC 4038/39 are 
listed in Table [2] and the slit positions are illustrated in 
Fig. [T] We note that the IRS position of the nucleus 
of NGC 4038 does not exactly coincide with the radio 
nucleus but is slightly shifted toward the north to also 
include the nearby, active star forming region. 

The data were pre-processed by the Spitzer Science 
Center data reduction pipeline version 13.2 (Observer's 
Manual, 2004). To avoid uncertainties introduced by the 
flat-fielding in earlier versions of the automated pipeline 
processing, we started from the two-dimensional, unflat- 
fielded data products. These products are part of the 
"basic calibrated data" (BCD) package provided by the 
Spitzer Science Center. 

First, all images were cleaned of bad pixels using the 
IDL procedure IRSCLEAN with a 'Maskval=128'. The 
various 2-dimensional spectra from the same nod posi- 
tion were median-combined within the Spectral Mod- 
elling, Analysis , and Redu ction Tool (SMART) version 
5.5.1. (|Higdon et al.|[200^ . A median 'sky' spectrum 
was computed from both nods of the sky position (#9) 
and subtracted from each median nod image. The spec- 
tral extraction was also done within SMART, using full 
aperture extraction. The 'fluxcon' tables were disabled 
by setting them to 'si'. Instead, the spectra were flat- 
fielded and flux calibrated by multiplication with the rel- 
ative spectral response function (RSRF) using the IRS 
RSRF of the standard star ^-Dra for both SH and LH 
spectra. We have not corrected for periodic "detector 
fringing" in the spectra since these artifacts have no ef- 
fect on the analysis carried out in this paper. 



Fig. 1. — Illustration of the IRS SH slit positions for t he given 
observing date, over-plotted on the IRAC Sfira map (Wa ng et al.l 
120041) of the Antennae. North is toward the upper right. 

After these steps there was a noticeable and expected 
mismatch between the fluxes measured by the IRS SH 
and LH modules at the same wavelengths. This mis- 
match is most likely due to extended emission, picked 
up by the wider LH slit. For an ideal point source there 
would be very little mismatch while for uniform surface 
brightness SH would only see 21% of the flux measured 
by LH. In principle there are three possibilities to match 
the two spectra: (i) scale LH down to match SH, (ii) 
scale SH up to match LH, or (in) use a wavelength de- 
pendent scaling factor that affects both SH and LH. (i) 
is best if the spectrum of a compact source surrounded 
by extended emission is of interest, (ii) should be used 
to quantify the extended emission, and (Hi) is arguably 
the best method to account for both components but re- 
quires additional assumptions. Since we are most inter- 
ested in the spectral information from unresolved SSCs 
- which are fully covered by both SH and LH slits - we 
have chosen method (i). In any case, the quantitative re- 
sults derived in this paper are from the lines covered by 
the SH module and are not affected by the scaling. The 
scaling factors 'SF' applied to the LH fluxes are listed in 
Table [2] The resulting IRS 'hires' spectra from the two 
nuclei and the six infrared peaks are shown in Fig. [3l 
The total observing time was 4.4 hours. 

2.2. 'Lores ' Observations 

We have also made observations in the IRS low resolu- 
tion ('lores') mode at i? = 65 — 130 with the Short-Low 
(SL) [5.2 - U.5^iTn] and Long-Low (LL) [14.0 - 38.0/im] 
modules. The angular slit sizes are listed Table [T] The 
observations were made with both 'lores' modules in 
spectral mapping mode. The AORs were intentionally 
designed to make use of the facts that the sub-slits SLl 
and SL2 are adjacent and record data simultaneously, 
and that the orientation would have changed after six 
months (half a solar orbit). However, an unfortunate de- 
lay of six months in the scheduling doubled the spatial 
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TABLE 2 

Parameters of the Pointed Observations with Cluster Designations from the Literature 



Name 


KA (JiOUOj 


Dec (J2U0U) 


LH SF^ 


Wby5 


W99'^ 


WZU2 


305"^ 


NOO^ 


nucleus 4038 


12 01 53.01 


-18 52 02.7 


0.77 












nucleus 4039 


12 01 53.54 


-18 53 10.2 


0.63 












peak 1 


12 01 54.98 


-18 53 05.7 


0.65 


80 




3 


157 


2-1 


peak 2 


12 01 54.58 


-18 53 03.4 


0.56 


86(89/90) 


11/9/14/30/33/39 


5 


136 


2-6 


peak 3 


12 01 55.39 


-18 52 48.9 


0.50 


119/120(/117) 


19/17 


7 


176 


4-2 


peak 4 


12 01 50.42 


-18 52 12.6 


0.59 


405 


2 


25 




11-2 


peak 5 


12 01 54.75 


-18 52 51.1 


0.45 


115 






148 


3-5 


peak 6 


12 01 54.80 


-18 52 13.5 


0.56 


384/382/389 


18 


24 


154 


5-5 


reference sky 


12 01 49.00 


-18 52 50.0 















^ Scaling factor to multiply with the LH fluxes to match SH.'' WS95: optical imaging with the HST 
WF/PC l|Whitmore fc Schweizer 19951).'^ W99: optical imaging with HST WFPC2 llWhitmore et al .1119991) .'' WZ02: 
IWhitmoreX Zhang (2002)= BOS: near-IR imaging with Palomar WIRC IIBrandl et al.ll2005l) .^ NOO: 4 and 6 cm radio 
observations with the VLA (.Neff fc Ulvestadll2000i) . 




Fig. 2. — Illustration of the areal coverage of the IRS SL spectral 
map overlaid on the IRAC 8/xm map. The larger part on top is the 
initial map while the lower, smaller area indicates the "filler" part. 
North is toward to upper right. 



offset between SLl and SL2 instead of providing the same 
spatial coverage. The problem was detected and cor- 
rected by requesting another SLl map six months later, 
specifically designed to fill in the missing parts. The 
mapped areas are illustrated in Fig. [5] The total map- 
ping time of all observations combined was 7.1 hours. 

All data were processed with the standard pipeline ver- 
sion 13. Further reduction and analysis started from the 
two-dimensional, flat-fielded data products, which are 
part of the "basic calibrated data" package provided by 
the Spitzer Science Center. 

Most of the reduction and spectral analysi s of the large 
data c ube has been made with CUBISM ([Smith et al.l 
dOOTi). CUBISM is an IDL program developed by the 
SINGS Legacy team to combine slit spectra at various 
angles and positions and create spectral maps and local 
spectra. Artifacts in the resulting spectral maps were 
cleaned with a combination of a local noise estimate 
and median replacement, and the procedure IRSCLEAN, 
version 1.3. The latter is based on a multi- resolution 



analysis algorithm (jMurtagh et al.|[l995f) . where an im- 
age is subsequently smoothed to different scales, using 
a wavelet transform. IRSCLEAN is a contributed IRS 
software, which we modified slightly to use a more lo- 
calized noise estimate on a sub-slit basis. We have also 
used Gaussian rather than uniform noise estimates. The 
algorithm was first tested on the "sky" spectra. Fig. d] 
shows the resulting line maps for six important spectral 
features: the fine-structure lines of [Ne II] and [S IV] , the 
11.3/Ltm and 8.6/j,m emission feature of PAHs, and the 
molecular hydrogen H2 S(2) and S(3) lines. 

3. ANALYSIS 

Fig. [3] reveals the spectral richness of the high signal- 
to-noise spectra. Important common features include the 
PAH emission bands, silicate absorption features, ionic 
fine-structure lines, atomic and molecular hydrogen lines, 
and the spectral continuum. In this Section we describe 
how the quantities relevant to our discussion have been 
derived from the 'lores' and 'hires' spectra. We empha- 
size the complementary character of the 'lores' and 'hires' 
data in our analysis. While the latter is used to derive 
accurate line fluxes the former is mainly used for a qual- 
itative assessment of the spatial characteristics and for 
cross-checks with the 'hires' data. 

3.1. Continuum Flux Densities 

In order to characterize the slope of the spectral contin- 
uum we have derived the flux densities F^, at two narrow 
wavelength intervals, namely at 14.75 — 15.25/xm and at 
29.5 — 30.5/im, here referred to as i^is^m and Fao^im, re- 
spectively. These wavelengths have been chosen as they 
are not significantly affected by known emission or ab- 
sorption features. Both intervals contain about 20 reso- 
lution elements, and the continuum within that narrow 
range is assumed to have a linear slope. The flux densi- 
ties of each resolution element within the wavelength in- 
tervals have been averaged to reduce their sensitivity to 
noise spikes and narrow spectral features, resulting in the 
quasi- monochromatic flux densities Fis^m and Fsopm- 

This appro ach is similar to the one described in 
iBrandl et all (p006i ) and is arguably the best direct esti- 
mate of the spectral continuum. The flux measurements 
are listed in Table [H The continuum fluxes will be dis- 
cussed in detail in Section [131 



3.2. Poly cyclic Aromatic Hydrocarbons (PAHs) 
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[S IV] [Ne II] [Ne III] [S lll][Ar lll][Fe lll][0 IV] [S lll][Si lit (S IV] [Ne II] [Ne III] [S lll][Ar lll][Fe lll][0 IV] [S lll][Si II] 



I I I I I I I I 

H, S(2)H(6-7) H,S(1) H, S(0) 

I I 




12 15 20 25 30 35 10 

Rest frame wavelength [^m] 



12 15 20 25 30 35 

Rest frame wavelength [jam] 



Fig. 3. 
in Fig.[T] 



— IRS SH+LH spectra of the two nuclei of NGC4038 and NGC4039 and the six infrared peaks listed in Table [2] and illustrated 
Most of the detected spectral lines are labeled. 
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TABLE 3 
Continuum Fluxes 







^30 fj,m 




Position 


[mjy] 


[mJy] 




nuc 4038 


103 


649 


0.159 


nuc 4039 


39 


334 


0.117 


peak 1 


383 


1994 


0.192 


peak 2 


260 


1853 


0.140 


peak 3 


79 


787 


0.101 


peak 4 


32 


231 


0.137 


peak 5 


69 


699 


0.099 


peak 6 


35 


233 


0.150 



Note. — The flux densities listed here and 
in the subsequent tables are the values directly 
measured within the IRS SH slit aperture of 
4."7 X ll."3. The fluxes longward of 20/xm were 
measured with the larger IRS LH slit but have 
been scaled down by the scaling factors listed 
in Table |2] 

The spectra in Fig. [3] show that the spectral shape at 
shorter wavelengths is largely dominated by the strong 
emission features from PAHs. The strength and equiva- 
lent width of the PAH fe atures were measure d using the 
IDL program PAHFIT (jSmith et al.ll2007H l. PAHFIT 
decomposes the spectra into the individual contributions 
from starlight, thermal emission from dust, atomic and 
molecular emission lines, and the PAH emission bands. 
The latter, combined with extinction dominated by the 
silicate absorption bands around 9.7/im and 18/im, is 
the main application of the routine and both are treated 
meticulously. 

We have chosen to measure the below listed PAH fea- 
tures from the 'hires' spectra rather than the 'lores' maps 
for two reasons. First, due to the longer integration times 
the signal-to-noise is higher than at individual positions 
of the spectral map. Second, all other spectral diagnos- 
tics were derived from the 'hires' spectra as well, and an 
accurate comparison does benefit from using the same slit 
apertures. However, since PAHFIT has been designed to 
analyze IRS 'lores' spectra and does not accurately fit 
narrow emission lines we removed the emission lines from 
the 'hires' spectra for the purpose of the PAH measure- 
ments. This was done by linear interpolation between the 
fluxes directly short-wards and long- wards of the emis- 
sion lines. PAHFIT uses Drude profiles to fit 15 PAH 
features of different width in the IRS-SH range. 

The thermal dust continuum is fitted by a combination 
of modified black bodies at fixed temperatures of 35, 40, 
50, 90, 135, 200 and 300 Kelvin. PAHFIT can automat- 
ically correct for extinction estimated from the depth of 
the 9.7/im silicate absorption feature. However, we have 
turned off this automatic feature since the 'hires' spectra 
cover only part of the 9.7/im silicate absorption and we 
noticed that, under these conditions, extinction is being 
used as a free parameter to optimize fit results without 
strong physical motivation. 

The central wavelength was fixed and the full width at 
half maximum (FWHM) was allowed to vary by at most 
20%. Several adjacent, broad PAH features cannot be 
resolved and are blended into one broad complex, and 
we list the combined fluxes from both components. Such 
cases are the features at 11.23 and 11.33/im, and at 12.62 
and 12.69//m. The PAH features at 14.04, 15.9, 18.92 
and 33.1/im were not included in the analysis because 



of their intrinsic weakness. The measured PAH fluxes 
and equivalent widths are listed in Tabled! The listed 
uncertainties include the uncertainty in the absolute flux 
calibration, which we estimate to be in the order of 10%, 
and the fit error as given by PAHFIT. Both contributions 
have been added in quadrature. 

The PAH strength in the 'lores' spectra has been visu- 
alized and analyzed within CUBISM in a more qualita- 
tive way. We co-added the emission from the spectral ele- 
ments covered under the emission feature and subtracted 
the continuum flux. The latter was derived from a lin- 
ear interpolation between the continuum fluxes to both 
sides of the feature, which were each computed from the 
median of typically three spectral resolution elements. 

We produced spectral maps in the 6.2/im, 8.6/im and 
11.3/im PAH emission features. The former closely re- 
sembles the one at 8.6/im and is not shown in this paper. 
Since the 12.7/im PAH map is affected by the [Ne II] 
line (see Section [373]) we did not include it in the further 
analysis. The 11.3/im feature is the one most affected by 
silicate absorption. At low spectral resolution the 8.6/im 
map has signifl cant contamination from the broad 7.7/im 
PAH feature (iSmith et al.ll2007bl ). However, since both 
species are similar in nature this will not affect the qual- 
itative conclusions drawn from the spectral maps. The 
8.6/im and 11.3/ini PAH maps are shown in Fig. (4] 

3.3. Fine- structure and Hydrogen Emission Lines 

As can be seen in Fig. [3] the mid-IR wavelength range 
contains numerous strong fine-structure emission lines, 
most prominently from sulphur and neon. We measured 
the line fiuxes from the 'hires' spectra using the Gaussian 
line fitting tool in SMART and a linear baseline fit. The 
results are listed in Table [5l The quoted uncertainties 
include the uncertainty in the absolute fiux calibration, 
which we estimate to be of the order of 10%, and the 
fit errors provided by SMART. Both contributions have 
been added in quadrature. 

Our measurements of [Ne III] / [Ne II] are in the range 
of 0.08 - 0.75, and hence lower th an th e [Ne III] / [Ne II] 
- 0.8 derived by Kunzc e t'^ (IT991 with ISO-SWS 
spectroscopy for the overlap region. However, our "line- 
luminosity weighted" average for the peaks 1, 2, 3, 5, and 
6 together is 0.571q 'J^, still lower than but close to the 
ISO-SWS value of 0'8. 

Fig. [3] also shows the lowest pure rotational transitions 
of molecular hydrogen, namely the 0-0 S(0) 28.21/im, the 
0-0 S(l) 17.03/im, and the 0-0 S(2) 12.28/im lines. In ad- 
dition we find the atomic hydrogen lines of Humphreys-a 
H(7-6) 12.37/im and the HI (8-7) 19.06/im. The lat- 
ter, however, has only been detected in the nucleus of 
NGC 4039. The fine fiuxes are fisted in Table M The 
quoted uncertainties include the uncertainty in the ab- 
solute flux calibration, which we estimate to be in the 
order of 10% for the 'hires' spectra and 20% for the spec- 
tral maps, and the fit errors provided by SMART, again 
added in quadrature. 

The line strength in the 'lores' spectra has been vi- 
sualized and analyzed with CUBISM in the same way 
as the PAH features before. We concentrate the analysis 
on the strongest fine-structure and hydrogen lines within 
the range of the IRS-SL module to have similar spatial 
resolution for all maps (the slit width increases by a fac- 
tor of two long-ward of 15/im). Hence, we present the 



Spitzer Spectroscopy of the Antennae 



7 



TABLE 4 

PAH FEATURE STRENGTHS LONGWARD OF 10/xM FROM THE 'HIRES' SPECTRA 







TT a 


Fl2.0Mm 


a 


Fl2.7fim^ 


^ 13.5/J.m 


Fl4.2/jm 


a 


F 


16.5/im 


a 


Fl7.4fjm^ 


Fl7.9pm^ 


Position 


EWlo.TMm'^ 


EWii.apm'^ 


EWia.OMm'" 




EWl3.5Mm'' 


EWl4.2Mm'' 


EWie.SMm'" 


EWl7.4Mm'' 


EWl7.9Mm'' 


nuc 4038 


0.73 ± 0.40 


18.02 ± 1.40 


9.41 ± 


65 


14.24 ± 1.36 


3.89 ± 0.46 


1.20 ± 


24 


1 


27 ±0 


20 


6.80 ± 1.40 


0.46 ± 0.63 




0.05 ± 0.01 


1.31 ±0.15 


0.67 ±0 


08 


0.97 ±0.14 


0.24 ± 0.04 


0.07 ±0 


02 





06 ±0 


01 


0.29 ±0.04 


0.02 ± 0.01 


nuc 4039 


0.53 ± 0.09 


9.91 ± 0.69 


4.45 ± 


25 


7.17 ±0.59 


1.97 ± 0.18 


6.06 ±0 


10 





51 ±0 


08 


3.53 ± 0.36 


0.30 ± 0.10 




0.09 ± 0.02 


1.70 ± 0.19 


0.75 ±0 


09 


1.18 ± 0.16 


0.31 ± 0.04 


0.09 ±0 


02 





06 ±0 


01 


0.42 ± 0.06 


0.03 ± 0.01 


peak 1 




8.49 ±0.90 


9.83 ±0 


98 


17.70 ±1.78 


15.40 ± 1.54 


8.06 ±0 


81 


1 


04 ±0 


11 


1.56 ±0.16 


0.92 ± 0.09 






0.20 ±0.02 


0.21 ±0 


02 


0.34 ±0.03 


0.27 ±0.03 


0.13 ±0 


01 





02 ±0 


01 


0.02 ±0.01 


0.01 ± 0.01 


peak 2 


0.67 ±0.25 


14.41 ± 1.57 


9.54 ±0 


77 


17.54 ±1.98 


7.57 ±0.76 


3.10 ±0 


53 





29 ±0 


36 








0.02 ± 0.01 


0.47 ±0.02 


0.28 ±0 


04 


0.48 ± 0.08 


0.19 ± 0.03 


0.07 ±0 


01 





01 ±0 


01 






peak 3 


0.40 ±0.11 


13.67 ±0.69 


6.42 ± 


34 


11.09 ±0.95 


3.30 ±0.28 


0.90 ±0 


19 





90 ±0 


15 


3.39 ±0.95 






0.04 ±0.01 


1.24 ±0.14 


0.56 ±0 


06 


0.92 ±0.12 


0.26 ± 0.03 


0.07 ±0 


02 





05 ±0 


01 


0.18 ±0.04 




peak 4 


0.68 ± 0.09 


12.17 ±0.82 


4.02 ±0 


56 


8.12 ±1.07 


1.42 ± 0.29 


0.67 ±0 


13 





58 ±0 


07 


3.43 ± 0.40 






0.10 ±0.02 


2.06 ±0.23 


0.75 ±0 


09 


1.60 ±0.24 


0.28 ± 0.04 


0.13 ±0 


03 





09 ±0 


01 


0.49 ± 0.07 




peak 5 


0.36 ± 0.19 


15.57 ±0.80 


7.06 ± 


84 


14.25 ± 1.11 


4.07 ±0.28 


1.46 ±0 


18 





90 ±0 


13 


4.03 ±0.78 






0.04 ±0.01 


1.79 ±0.21 


0.75 ±0 


08 


1.40 ±0.18 


0.37 ±0.05 


0.13 ±0 


02 





06 ±0 


01 


0.27 ±0.05 




peak 6 


0.70 ±0.08 


12.97 ±0.85 


4.18 ±0 


58 


9.00 ±1.17 


1.66 ± 0.30 


0.96 ±0 


11 





79 ±0 


09 


3.67 ±0.77 






0.10 ±0.02 


2.00 ±0.23 


0.70 ± 


08 


1.56 ±0.23 


0.29 ± 0.04 


0.16 ±0 


02 





11±0 


02 


0.50 ±0.05 





Note. — The listed measurements have not been corrected for extinction. Sec comment to Table[3]for sht aperture sizes. 
^ Integrated flux in units of lO'^^W cm Equivalent widtii in units of fiiri 



TABLE 5 

Fine-structure line fluxes in units of 10~^''Wcm~^ 



Line 


[S IV] 


[Ne II] 


[Ne III] 


[S III] 


[Ar III] 


[Fe III] 


[O IV] 


[Fe II] 


[S III] 


[Si II] 


Wavelength 


10.5lAtm 


12.8lAtm 


15.56/xm 


18.7lAtm 


21.83Atm 


22.93Aim 


25.89Aim 


25.99Aim 


33.48AJm 


34.82^im 


EP=- 


34.97 eV 


21.56 eV 


40.96 eV 


23.34 eV 


27.63 eV 


16.19 eV 


54.93 eV 


7.90 eV 


23.34 eV 


8.15 eV 


nuc 4038 




3.83 ± 0.40 


0.29 ±0.03 


1.89 ±0.20 




0.12 ± 0.02 


0.12 ± 0.04 


0.21 ±0.03 


2.75 ±0.28 


3.88 ±0.39 


nuc 4039 


0.34 ±0.04 


2.41 ±0.25 


1.14 ±0.12 


1.91 ±0.20 




0.11 ± 0.01 


0.10 ±0.69 


0.23 ± 0.31 


2.67 ±0.27 


3.25 ±0.33 


peak 1 


1.25 ±0.14 


7.05 ±0.77 


4.44 ±0.46 


5.14 ±0.52 


0.05 ±0.08 


0.22 ± 0.04 


0.14 ±0.03 


0.16 ± 0.04 


7.64 ± 0.80 


4.73 ±0.51 


peak 2 


1.87 ±0.19 


6.64 ±0.70 


4.86 ±0.49 


6.22 ±0.69 




0.30 ±0.09 




0.26 ± 0.14 


8.72 ±0.93 


5.19 ±0.54 


peak 3 


0.78 ±0.08 


3.98 ±0.41 


2.48 ±0.26 


3.21 ±0.32 




0.18 ± 0.02 


0.13 ± 0.02 


0.21 ± 0.03 


6.09 ±0.62 


4.90 ±0.50 


peak 4 


0.05 ±0.01 


2.11 ±0.25 


0.65 ±0.07 


1.37 ±0.14 




0.06 ± 0.01 


0.05 ± 0.01 


0.10 ±0.01 


3.18 ±0.32 


3.18 ±0.32 


peak 5 


0.39 ±0.04 


3.33 ±0.34 


2.31 ±0.23 


2.43 ±0.25 


0.05 ±0.02 


0.17 ±0.02 


0.17 ±0.04 


0.22 ± 0.02 


5.14 ±0.53 


4.51 ±0.46 


peak 6 


0.14 ±0.02 


2.47 ±0.27 


0.95 ±0.10 


1.71 ±0.18 


0.02 ±0.00 


0.07 ±0.01 


0.06 ± 0.01 


0.10 ±0.01 


2.55 ±0.26 


3.23 ±0.33 



Note. — Tlic listed measurements have not been corrected for extinction. See comment to Table [3] for slit aperture sizes. 
^ Excitation potential to create this ion. 



TABLE 6 

Line fluxes of molecular and atomic hydrogen in units of 10~^°Wcm~^ observed in the IRS SL 

AND SH modules. 



Transition 


H2S(3)'' 


H2S(2)'' 


H2S(2)'^ 


HaSCl)*^ 


H2S(0)'^ 


H (7-6)'' 


H (8-7)'' 


Wavelength 


9.67/im 


12.28/im 


12.28AJm 


17.03Ajm 


28.22fjm 


12.37Ajm 






1015 K 


510 K 


510 K 


170 K 


K 






nuc 4038 


0.71 ±0.15 


0.52 ± 0.11 


0.64 ± 0.08 


1.43 ± 0.15 


0.72 ±0.21 


0.06 ±0.01 




nuc 4039 


3.76 ± 0.76 


1.99 ± 0.42 


1.11 ±0.12 


2.37 ±0.24 


0.43 ± 0.08 


0.04 ±0.05 


0.05 ±0.01 


peak 1 


< 0.50<* 


0.39 ± 0.09 


0.62 ± 0.08 


1.64 ±0.17 


0.81 ±0.09 


0.15 ±0.02 




peak 2 


0.43 ±0.11 


0.45 ± 0.11 


0.36 ± 0.06 


1.23 ±0.13 


1.09 ±0.24 


0.11 ±0.04 




peak 3 


0.53 ±0.12 


0.45 ± 0.10 


0.66 ± 0.07 


1.77 ±0.18 


0.90 ±0.10 


0.11 ±0.03 




peak 4 






0.11 ± 0.03 


0.19 ± 0.02 


0.22 ±0.05 


0.06 ±0.02 




peak 5 


0.69 ±0.14 


0.67 ±0.15 


0.95 ± 0.11 


2.38 ± 0.24 


0.87 ±0.10 


0.11 ±0.02 




peak 6 






0.13 ± 0.03 


0.24 ±0.03 


0.16 ±0.02 


0.05 ±0.01 





Note. — The listed mea sure ments have not been corrected for extinction. See comment to Table |3] for the 
'hires' slit sizes and Section 13.31 for the 'lores' apertures. 

^ Observed with the IRS-SL module. The quoted uncertainties only include the fitting and photometric 
uncertainties but no systematic errors of the spectral mapping.'' Observed with the JRS-SH module.'^ Energy 
of lower level in temperature equivalents.'' Upper limit only, due to an artifact in the spectral map. 




Fig. 4. — IRS-SL spectral maps at the [Ne II] and [S IV] lines (top), the ll.S/^m and 8.6/xm PAH features (center), and the S(2) and S{3) 
rotational lines of H2. The spectral continuum has been subtracted and the images shown represent the pure line intensities. The maps 
are shown in hneax scaline. North is ud. east is left. The imaees are aoDroximatelv 135" (lAA koc) aside. For a discussion see Section 14.41 
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[S IV] map in Fig. [4] instead of [Ne III] which hes outside 
the range covered by the IRS-SL module. 

We note that in 'lores' mode, the strong, adjacent emis- 
sion features of [Ne II] at 12.81/im and PAH at 12.7/im 
cannot be fully spectrally resolved. From the 'hires' spec- 
tra we know that, on average, the luminosity in the PAH 
feature is about three times higher and thus a serious po- 
tential contaminant. A cross-contamination may lead to 
an over-estimation of the line flux, but also to an over- 
estimation of the continuum flux. For that reason we 
have excluded the 12.7/im PAH feature from our analy- 
sis. A comparison of the [Ne II] line flux in the nucleus of 
NGC4039 with the 'hires' spectrum, where the compo- 
nents are clearly resolved, indicates that the PAH con- 
tamination to the line flux is small (~ 20%), and that 
the [Ne II] map is suitable for qualitative arguments. 

The 0-0 S(l)17.03//m, the 0-0 S(2)12.28Afm and the 0- 
S(3)9.66/im emission lines of molecular hydrogen were 
also detected in the 'lores' spectral maps. The S(2) and 
S(3) lines were both measured with the IRS-SLl mod- 
ule, while the S(l) line was measured with the IRS-LL2 
module at ^ 2x lower spatial resolution. Unlike the 
fixed 'hires' slit sizes and orientation, the apertures for 
the spectral maps, from which the 'lores' fluxes were ex- 
tracted, were best matched to the size of the emitting 
region and vary from 5x4 pixels (960 x 768 pc) to 10x11 
pixels (1.92 X 2.11 kpc) for the area around the nucleus 
of NGC4039. The fluxes of the S(2) and S(3) lines have 
been added to Table \6\ Fig. [4] shows the spectral maps 
at the [Ne II], [S IV], and the H2 0-0 S(2)12.28/im and 
S(3)9.66/im lines. 

3.4. Deriving H2 Temperatures and Masses 

The H2 lines serve as crucial diagnostics of the condi- 
tions of the ISM. The H2 temperatures were calculated 
independently for the 'lores' and the 'hires' data as de- 
scribed in more detail in the Appendix and are listed 
in Table [9l Values range from approximately 270 K to 
370 K. 

However Table [9] does not list the uncertainties in the 
temperature and mass estimates because they are domi- 
nated by three components, which we discuss here. First 
are the errors in the line flux measurements, listed in 
Table [H These are typically about 3%. Second, the 
different clusters suffer different amounts of extinction. 
Correcting for extinction has two effects: it increases the 
line fluxes and it changes the line ratios. For instance, 
a "typical" optical depth of T9.7 — 0.2 would raise the 
temperature determined from the S(2) and S(3) lines by 
about 4% (~ 12 K) since the S(3) is located in the center 
of the silicate band and suffers most from extinction. If 
we use the S(l) and S(2) lines instead, the temperatures 
rises only by 0.7% 2K). The effect of extinction on 
the mass estimate is more complicated since it increases 
linearly with the corrected line flux but depends also on 
the temperature. For rg.7 = 0.2 the H2 mass increases by 
approximately 5%. Given the large uncertainties in the 
Av or T9.7 for the clusters, as discussed in Section 13.51 
and summarized in Table [7l we decided not to correct 
our tabulated estimates for extinction. Third, aperture 
effects play an important role. The 'hires' aperture is 
defined by the slit width and length and not matched to 
the size of the H2 emitting region. For the 'lores' values 
the aperture was matched to the size of the region (Sec- 



TABLE 7 
Extinction estimates 



Position 






Av 




peak 1 


4.23 


6.2 ±0.3 


4.3 ±0.3'= 


0.19 


peak 2 


0.18 


0.7 ±0.1 


1.4 ±0.3'= 


0.13 


peak 3 


0.14 


4.8 ±0.4 






peak 4 






0.6 ± 0.3'= 




peak 5 


11.81 


10.3 ±0.5 




1.03 


peak 6 






0.72^^ 


0.05 



^ Av from a comparison of ncar-IR photometry of 
various broad- and narrow^-band fluxes with STAR- 
BURST99 l IMengel c t ai;'200 ^> Ay from Br-^/Pa g and 
the Brackett series (Sniidcrs van der 
from Br^/H„ I M cngel et"an 120011) .'^ from H^/Hfj 

l^astian ct al. 200^.^ rg y estimates from baseline fitting 
f5cction l3.5t . 

tion l3.3|) . This systematic error is refiected in Table [9] by 
the difference between the corresponding mass estimates, 
and is the dominating uncertainty. 

We have also calculated the H2 masses following the 
procedure outlined in the Appendix. The derived masses 
are listed independently for the 'lores' and the 'hires' 
data in Table [51 Again, the two independent mass esti- 
mates using different lines and measurements from differ- 
ent instruments agree within 15 % (except for peak 3), 
which provides confidence in the reliability of our mea- 
surements. A discussion of the derived H2 masses is given 
in Section |4TT1 

3.5. Extinction Estimates 

The wavelength range covered by the ffiS-SL and SH 
modules includes two important diagnostics to estimate 
the amount of dust: a broad Si=0 stretching resonance, 
peaking at 9.7/im, and an even broader 0-Si-O bending 
mode resonance, peaking at 18.5/im. We have estimated 
the extinction using the apparent optical depth in the 
9.7/im silicate feature. We assume extinction from a fore- 
ground screen, which is certainly an over-simplification. 
However, unlike for starburst galaxies as a whole, it is a 
reasonable approximation for the extinction toward indi- 
vidual, embedded SSCs, where the luminosity is provided 
by a central stellar cluster which is surrounded by dust 
and gas clouds and PDRs. 

The depth of the 9.7/im silicate feature was directly 
measured as the natural logarithm of the ratio of ob- 
served fiux to the nominal mid-IR continuum at 9.7/im. 
The latter can be best estimated (for PAH dominated 
spectra) by a power law fit anchored at 5.5/im and 
14.5/im. The metho d is described and illustrated in 
Fig. 2 of ISpoon et al.l (j2007). Our estimates of the opti- 
cal depth are listed in the right column of Table [T] 

We disabled the PAHFIT default setting to automati- 
cally correct the measurements for extinction, and hence 
the values in Table [4] are the uncorrected, directly mea- 
sured fiuxes. This has been done for three reasons. First, 
PAHFIT has many free parameters, and given the lim- 
ited spectral range of the JRS-SH module, which does 
not fully cover the 9.7/im silicate band, extinction could 
be used as yet another free parameter to optimize the 
overall fit. Second, we do not, a priori, know whether a 
dust screen or mixed geometry provides a more realistic 
description. Third, the extinction at mid-IR wavelengths 
is generally rather low in these objects. 
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Brandl et al. 



For a large sample of 172 near-IR sources in the An- 
tennae, iBrandl et all (|2005| ) found that the average ex- 
tinction is about Ay ~ 2, while the reddest clusters may 
be attenuated by up to 10 magnitudes. From ISO-SWS 
spectroscopy of the hydrogen recombination lines Br-y 
and Br^ in the overlap region iKunze et alJ (|1996f l de- 
rived an extinction of Ay ~ 70 (mixed case) or Ay ~ 15 
(screen extinction). The literature values are compared 
with our measurements in Table [T] 

The conversion from T9.7 to Ay strongly depends on 
the assumed extinction law. Toward the G alactic Center 
there are various estima tes: Ay/rc, 7 ~ 6.7 (jMoneti et al.1 
[200l. Ay /to 7 ~ 7.8 (lL^ [l99m. and Ay/rc, 7 ~ 9.0 
(iMathisI [199a ). For the Solar neighbourho od iMathisI 
(|l990^ derived Ay/rg.y - 18.5, and iDr^^i (|1989f ) cal- 
culated Ay/rg.y ^ 18.1. Given this broad range of con- 
version factors, we find that the estimates for T9.7 from 
baseline fitting agree reasonably well with the other val- 
ues listed in Table [71 

We note that the Tg,7 values indicate less extinction 
than derived from observations at shorter wavelengths, 
and the ratio of Ay /tq,7 varies significantly from source 
to source. This may be caused by a dilution effect due to 
the large area covered by the IRS slits. An extended 
component of strong mid-infrared continuum emission 
will decrease the relative depth of the silicate absorp- 
tion feature. This effect is likely most predominant for 
peak 1. Furthermore, optical methods mainly probe the 
absorption by graphite-based dust particles, while the 
mid-IR methods are mainly sensitive to distinct silicate 
absorption. In the diffuse ISM graphite and silicate- 
based dust are uniforin ly mixed an d tightly cor r elated 
(jRoche fc AitkenI I1984D . Recently, iChiar et all (|2007f ) 
found that in dense clouds with Ay > 12 the linear rela- 
tion between optical and mid-IR estimates breaks down 
and T9.7 underestimates the real amount of dust. How- 
ever, most of our regions are below Ay ^ 12 and the 
discrepancy between the extinction values in Table [7] is 
primarily given by systematic uncertainties rather than 
ISM chemistry. 

4. RESULTS AND DISCUSSION 

In this Section we discuss and interpret our observa- 
tional findings from both the spectral maps and the spec- 
tra of the two nuclei and six infrared peaks. Our aim is to 
characterize the conditions under which SSCs form and 
evolve, their properties, and how their presence affects 
the surrounding interstellar medium. We start with a 
discussion of the individual regions, discuss qualitatively 
the ISM properties across the central region of the An- 
tennae, investigate dust temperatures and radiation field, 
and derive cluster masses and star formation rates. We 
discuss the strength and variability of the PAH features 
and focus specifically on the strength of the H2 emission 
and the temperature and excitation of the molecular hy- 
drogen in the Antennae. 

4.1. The Nuclei of NGC4038 and NGC4039 

Bas ed on sub-arcsecond near-IR spectroscopy from 
Keck, [Gilbert et al.l (12000*) found that the spectrum of 
the nucleus of NGC4039 is marked by a strong stellar 
continuum and bright, extended H2 emission. These au- 
thors also found strong photosphcric absorption of Mg I, 
Na I, and Ca I as well as from the CO Ah' = 2 band 



head absorption, indicating that the continuum is domi- 
nated by old giants and red super-giants. The Br^ line, 
an indica tor of recent massiv e star formation, was not 
detected. iMengel et al.l ()2001[ ) found starburst activity a 
couple of arcseconds north of the nucleus of NGC4038, 
con sistent with an age o f around 6 Myr. For the two nu- 
clei iMengeTeFal] ()2001[ ) estimated an age of 65± 15 Myr 
from CO absorption. It is clear that two galactic nu- 
clei ha ve lower star form ation activity than the overlap 
region (jWang et al.ll2004[ ). 

Chandra detected several X-ray sources near the ra- 
dio positions of the nuclei of NGC 4038 and NGC 4039 
(jZezas et al.ll2002aD . The X -ray sources identified with 
the nuclei are both luminous and spatially extended. The 
NGC 4038 nucleus spectrum is very soft and steep, and 
consistent with thermal emission by a supernova-driven 
super- wind, possibly with some contribution by X-ray bi- 
naries (jZezas et al.H'2002allbf ). The nucleus of NGC 4039 is 
more luminous in X-rays, but its spectrum also suggests 
a combinati on of X -ray binaries and compact supernova 
remnants (jZezas et al. 2002a bl. The N GC 4039 nucleus 
also has a very steep radio spectrum ()Nefl[ fc UlvestadI 
|2000[) suggesting that its radio emission is non-thermal, 
arising from supernova remnants. These observations are 
consistent with the intermediate starburst age of about 
65 Myr estimated by Mcngcl ct al. (2001). The super- 
novae may be responsible for heating the ISM in the nu- 
clear region to the observed X-ray temperatures, and the 
X-ray emission, and shocks arising from the supernovae 
may also be responsible for the high-excitation line emis- 
sion. At any rate, neither X-ray nor radio observations 
have provided strong evidence for activity due to black 
hole accretion in either nucleus. 

The IRS spectra of the two nuclei are quite similar 
in terms of their PAH strength and silicate absorption 
(Fig. [3]). However they differ significantly regarding their 
mid-IR fiuxes, 20-30/im slopes, and fine-structure lines. 
As mentioned in Section [^TT] the observed position of the 
nucleus of NGC 4038 is slightly offset to the north of 
the radio nucleus, and the wide IRS-SH slit includes the 
nearby star forming region as well. Hence, the detected 
mid-IR emission is likely affected by the starburst ac- 
tivity to the north of the nucleus of NGC 4038. The 
nucleus and the nearby starburst region can be seen in 
most of the spectral maps in Fig. 3] as an extended struc- 
ture with two separate components. The strong high 
excitation lines in the NGC 4039 nucleus pose somewhat 
more of a puzzle. The [Ne III] / [Ne II] ratio is almost 
six times higher in NGC 4039, and the [S IV] line is also 
very strong in NGC 4039, while it has not been detected 
in NGC 4038. 

In the mid-IR AGN activity can usually be traced 
by the h igh excitation [O IV] and [Ne V] emis- 
sion fines (fSpinoglio fc Malkanlfl99HSturm et "all 120021 : 
iWeedman et al.|[2005[ ). However, the best indicators, the 
[Ne V] lines at 14.3 and 24.3/im, were not detected in 
either nucleus. The [O IV] line is not produced in mea- 
surable quantities in Hii regions but can be excited by 
hot Wolf-Rayet stars, sho cks, and AGN. In a study of 
starburst galaxies with ISQ iLutz et al.l (|1998[ ) defined the 
ratio [O IV] / ([Ne II]-h0.4x[Ne III]), and found values 
ranging between 0.006 (M82) and 0.062 (IIZw40) with a 
mean of 0.022. As can be seen from Table[5]our ratios for 
the peaks 1, 3, 4, 5 and 6 are 0.016, 0.026, 0.021, 0.040 
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- NGC 4039 nucleus 
■ Starburst model 




20 25 
Rest frame wavelength [/im] 

Fig. 5.— IRS spectrum of the NGC4039 nucleus (thin 
curve) with the be st fitting starburst SED model template from 
IGroves et al.l II2008I ) (thick dashed curve). 

and 0.021, respectively - all in excellent agreement with 
the mean starburst value of iLutz et alj ()1998D . For the 
two nuclei of NGC 4038 and NGC 4039 we find ratios of 
0.030 and 0.035, respectively. Both ratios are also well 
within in the range covered by "pure" starburst systems 
and provide no evidence for AGN activity. This find- 
ing is further supported by the large equivalent widths 
of the PAH features (Tab. [3]), which would not be ex- 
pecte d in the surrounding s of an accreting black hole 
(e.g.. lWeedman et al.|[2005[) . 

Besides the high excitation lines, arguably the best 
mid-IR diagnostic for the presence of an AGN is the 
spectral continuum. Hot dust around an AGN will sig- 
nificantly flat t en th e continuum slope (see e.g., Fig. 8 in 
IBrandl et all (f2006l) ). In Fig. [5] we compare the spec- 
tr um of the nuc l eus o f NGC 4039 to a starburst model 
of IGroves et al.l (|2008D chosen to best match the dust 
continuum and the PAH features. The determined pa- 
rameters of the fit indicate that the nuclear abundances 
are similar to Solar IZq), with the star formation 
occurring in somewhat distributed manner (logC = 4.5) 
in regions with standard ISM pressures Po/k = 10^ 
Kcm~^), but still reasonably embedded within their 
molecular gas birth clouds .fp DR ~ 1-0). (For a fu ll 
description of the parameters see IGroves et al.l ()2008D ). 
The good match between the observed continuum and 
a "pure" starburst model shows that the spectrum of 
the nucleus of NGC 4039 is fully consistent with emis- 
sion from dust that is heated by star formation only. 

4.2. The Properties of Individual SSCs 

Four of our six infrared peaks, namely peaks 1, 2, 
3, and 5, are located in the overlap region (Fig. [T|), 
which is the re gion of the most intense star formation i n 
the Antennae (jMirabel et al.lll998Hvrgroux et al.lll996D : 
peaks 4 and 6 are located outside the overlap region. In 
this subsection we summarize briefly some of their prop- 
erties based on information from the literature. Some 
cluster masses are discussed in Section [4.61 

Peak 1 coincides with a clump of molecular gas 
as traced by CO (peak CO-S and SGMC 4-5, 
IStanford et all [l990l : IWilson et al.] [2000l l200l . It is 
also the brightest source at centimeter wavelengths, 
and the spectral index of the radio emission indi- 
cates that the spectral energy distribution is dom- 



inated by the thermal radiat ion from Hii regions 
(jHummel fc van der HulstI I1986D with an equivalent of 
~ 5000 05 stars (iNeff fc UlvestadI [2000l) . However, 
the optical counterpart to this source corresponds to a 
very inconspicuous, faint, red source in the HST images 
(jWhitmore fc Schweizeil Il995l source 80) obscured by 
large amounts of gas and dust. Based on sub-arcsecond 
near- and i nid-IR spectro s copic data, both lGilbert et all 
(|2000( l and lSnijders et al] (|2007t l find an ionized gas den- 
sity of 10^ cm~'^. 

Peak 2 is the mid-IR counterpart to a bright and blue 
complex of clusters, whic h contains eight optical source s 
within a region of 1'.'5 (|Whitmore fc Schweized [l995h . 
corresponding to 160 pc in projection at our adopted dis- 
tance. Peak 2 also corresponds to the second brightest 
radio source in the Antennae. Like peak 1 it has a shallow 
spectral slope, typical for t he emission from Hil regions 
(|Hummel fc van der Huls^ [1981 iNeff fc UlvestadI [20Ml 
with an equivalent of ~ 3000 05 stars. The density of 
the ionized gas has been estimated to be ~ 10^ cm~'^ 
(jSniiders et al.ll2007[ ). 

Peak 3 corresponds to the center of the largest concen- 
tration of molecular gas in the overlap regi on. The coinci- 
dent sources CO J and part o f SGMC 1 (IStanford et all 
119901: IWilson et al.ll2000l. l2003l) contain a total molecular 
gas mass over 6-10® Mq . The slope of the correspond ing 
bright, compact radio source (|Neff fc UlvestadI l2000l ) is 
steep, arguing for a somewhat older cluster with its radio 
emission being dominated by supernova remnants. 

Peak 4 is located to the far west of the nucleus 
of NGC 4038. The radio observations show a mod- 
erately bright, compact source with a steep spec- 
tral slope, indicative of an evolved stellar populatio n 
(jHummel fc van der HulsillTMl INeff fc Ulvestadl[200l . 

Peak 5 is located in the overlap region, approxi- 
mately 1.6 kpc to the north-west of peak 1, and is 
the brightest sub-millimeter peak. It coincides with the 
sources CO.W and part o f SGMC 1 (|Stanford et al.lll990l : 
IWilson et al.][2000l ^OOl , which are huge concentrations 
of approximately 3.9 • 10* Mq of molecular gas. With 
an Ay ^ 11.81 it is one of t he most heavily em bed- 
ded clusters in the Antennae (jMengel et al.l [20051 ). Its 
op tical counterpart is identified as the red source 115 
in iWhitmore fc Schweized ()1995( ) . Pea k 5 also coincides 
with a bright, compact radio source (|Neff fc UlvestadI 
l2000f ). The spectral index suggests that the dominant 
radio emitters are supernova remnants. 

Peak 6, finally, is a mid-IR source located in a region 
of more extended, fuzzy IR emission, to the east of the 
nucleus of NGC 4038. The spectral slope of the radio 
emission is at the division point between the radiation 
field of Hii regions and supernova remnants. 

4.3. Cluster Ages 

The age estimates depend significantly on the method 
being used, and show significant scatter. The literature 
values on cluster ages of all our clusters are summarized 
in Table [SI The various age diagnostics indicate that the 
luminous IR sources in the northern part of the overlap 
region are somewhat older than the ones in the more 
active, southern overlap region. 

In Fig. [Sjwe compare the average cluster ages from the 
literature (Tab. [5]) to the radiation hardness measured 
by the [Ne III] / [Ne II] line ratio. There is a clear trend 
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TABLE 8 
Cluster Ages from the Literature 



Age'' 


Age" 


Agc'^ 


Aged 


Age<= 


Position [Myr] 


[Myr] 


[Myr] 


[Myr] 


[Myr] 



peak 1 


2.0 


2.3 


- 4.0 


3.5 


< 2.5 


2.5 


peak 2 


3.8 


2.3 


- 4.0 




< 3 


3.0 


peak 3 


2.0 


3.2 


- 4.9 


3.9 


< 3 


3.0 


peak 4 


7.0 




> 6 






7.0 


peak 5 




3.7 


- 5.1 


5.7 


3-5 


4.5 


peak 6 


8.4 


4.3 


- 5.7 






6.0 



^ from H„ llWhitmore fc Zhand 120021) '' from Br., equiva- 
lent width JMe ngel et aI.| |2005D ^ from Br., equivalent widtii 
JGilbert Gra liam 20071)*^ from the equivalent width of sev- 
eral hydrogen recombination lines JSniiders et al.ll2007D ^ av- 
erage age from the listed estimates 

of older clusters displaying a softer radiation field due to 
the disappearance of the early-type O and Wolf-Rayet 
stars after about 5 IVIyr. While this result may not come 
as a surprise it is reassuring that completely independent 
methods are consistent. 
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Fig. 6. — The hardness of the radiation field (as measured by 
the fine flux ratio of [Ne III] / [Ne II]) versus the average age of the 
SSCs from Table \E\ The horizontal dotted lines indicate the age 
ranges covered by the various estimates. 



We note that starburst activity exists to both east (in- 
cluding peak 6) and west (including peak 4) of the nu- 
cleus of NGC4038 at almost the same level, in terms 
of the [Ne II] and PAH intensity. In fact, peaks 4 and 
6 have almost identical mid-IR spectra (Fig. [3]), yet are 
separated by more than 6.6 kpc in projection. With a ve- 
locity dispersion of ~ 10 kms~^ (jWhitmore et al.ll2005f ) 
the dynamical time to connect the two points would be 
647 Myr, about two orders of magnitude longer than the 
cluster ages. This suggests that the dominant mode of 
triggering the formation of SSCs does not propagate hy- 
drodynamically but is governed by local processes. 

4.4. Physical Conditions across the Antennae 

In this subsection we discuss the spatial variations in 
the ISIVI as revealed by the spectral maps shown in Fig. [H 
These maps provide a more global picture of the physi- 
cal conditions across the central region of the Antennae 
galaxies. Unfortunately, the [Ne III] line, which we uti- 
lize heavily in Section B31 is not covered by the IRS-SL 
module and would have to be observed through the much 



wider IRS-LiL slit. However, given that the excitation po- 
tential of [Ne III] is reasonably close to that of [S IV] (see 
Table [5]) , we substitute the [S IV] line as a tracer of the 
harder radiation. 

Fig. m (top) compares the maps of [Ne II] with [S IV] . 
It takes only 21.6 eV to excite [Ne II], and the emis- 
sion traces primarily the regions where star formation 
occurred in the past ~ 10 Myr. (As a rule of thumb 
we assume that [Ne II] is primarily associated with stars 
younger than ~ 10 Myr, while [S IV] traces stars younger 
than ~ 4 — 6 Myr in a starburst). The [S IV] emission 
originates predominantly from two clusters peaks 1 and 
2 in the overlap region, with peak 2 being even brighter 
in [S IV] than peak 1 (Fig. 3]) . Although this has already 
been indicated by the spectra in Fig.[3l and qualitatively 
known since ISO (Mirabel e t al. 1998), it is once again 
remarkable how the bulk of the current massive star for- 
mation in an interacting system, which extends over tens 
of kilo-parsecs, is confined to just two compact regions. 
We also detect some [S IV] emission from the nucleus of 
NGC4039, but none from 4038 (Section [HI]). The [Ne II] 
map is also dominated by the two clusters peak 1 and 2, 
although to a much lesser degree than the [S IV] map, 
and also shows emission from the nucleus of NGC4038. 
Noticeable [Ne II] emission is also detected to the west 
of the nucleus of NGC 403 8, where numerous older SSCs 
reside (|Mengel et al.ll2005D . 

Overall, the PAH maps in Fig. 2] (center row) resemble 
the [Ne II] map^^ , but show noticeably more diffuse emis- 
sion. In "normal" galaxies a predominant fraction of the 
total PAH emission arises from the diffuse ISM, heated 
only by the interstellar rad iation field, with a only a low 
fraction rising from PDRs praine et al.ll2007f) . However 
in starbursting galaxies like the Antennae the situation 
is reversed, with the major part of the total emission 
coming from the regions of recent star-formation. Here 
the PAH maps trace primarily the PDRs and hence the 
environment of OB clusters. However, the fact that the 
[Ne II] and PAH map resemble each other qualitatively 
implies that we do not resolve the Hii region/PDR in- 
terfaces at the spatial resolution of the IRS-SL spectral 
map, which is between 200 - 380 pc, depending on wave- 
length and orientation. 

The most obvious difference between the [Ne II] and 
the PAH maps is the strength of peaks 1 and 2 relative 
to the nucleus of NGC 4038. The latter is the dominant 
source of PAH emission, while the embedded SSCs play 
only a minor role. Within the uncertainties, which are 
dominated by the lower signal-to-noise in the 6.2^m map, 
we find no significant differences between the 11.3/im and 
the 6.2^m PAH maps. (We note that the PSF in the 
former is also 1.8 times larger). The similarity between 
the two PAH maps agrees with our study of the PAH 
spectrum in Section 

The bottom row in Fig. [4] shows the spectral maps 
in the H2 S(2) and S(3) lines. Both maps are very 
similar. We find that approximately 45% of the total 
H2S(3) emission comes from a slightly extended region 
around the southern nucleus, NGC 4039. The nucleus 
of NGC 4038 does show some H2 emission, but very lit- 
tle compared to NGC 4039. The emission peaks of the 

^■^ See section 13.31 for a discussion of the contamination of the 
[Ne II] map by 12.7/im PAH emission. 
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clusters in the overlap region are more compact and also 
weaker. Most interesting here is the small spatial off- 
set of about two pixels, corresponding to approximately 
380 pc, that exists between the peaks of H2 emission and 
the spectral continuum, for both peak 5 and the nucleus 
of NGC 4039. The H2 S(2) emission peaks to the south in 
both objects and appears also, at least near the nucleus 
on NGC 4039, more extended. Since the offsets are with 
respect to the continuum emission, which was simulta- 
neously observed with the line emission, they cannot be 
due to spacecraft pointing errors or other observational 
effects. Rather they must arise from a physical displace- 
ment between the massive SSCs and the source of H2 
emission. 

Finally, we compared t he 8.6Atn i PAH map with the 
IRAC 8^im image of Wang et al.l ()2004[ ). Despite the 
slightly better spatial resolution of the IRAC image both 
maps look very similar with one exception: the emission 
from peak 1 is noticeably brighter in the IRAC image 
than in the 8.6/im PAH map. The spectra in Fig. [3] 
and the values in Table 2] show that the PAH equiva- 
lent widths of peaks 1 and 2 are significantly below those 
of the other clusters. Although the IRAC 8/im band is 
generally dominated by PAH emission, a proper contin- 
uum subtraction is obviously important for quantitative 
results. 

4.5. Dust Temperature and Radiation Field 

The mid-infrared continuum > 20/im is produced by 
the thermal emission from dust heated by the stellar clus- 
ters. The slope of the continuum reflects the temperature 
distribution of the dust. The presence of a hot dust com- 
ponent (100 K < T < 250 K) flattens the continuum slope 
in the wavelength range of 15 — 35/im, while a cold dust 
component (T < 50 K) steepens the slope. The contin- 
uum fluxes at 15/im and 30/im and their ratios are listed 
in TableO Peaks 1 and 2 have the shallowest slopes, indi- 
cating the highest dust temperatures, while peaks 3 and 
5 show the steepest continua, evidence for a dominating 
dust component at lower temperature. The difference is 
also illustrated in Fig. [7] for the most extreme cases of 
peaks 1 and 5, in com parison to th e aver age star burst 
galaxy spectrum from IBrandl et al.l (j2006f ) . The slope 
of the starburst template lies between the two Antennae 
clusters. 

Generally, the dust temperature depends on a com- 
plex relation between ISM density, grain size distribu- 
tion, stellar properties, strength of the radiation field, 
and the local geometry of cluster, Hii region, PDR, and 
surrounding ISM. We have investigated possible corre- 
lations of the dust temperature with cluster properties 
from the literature and our measurements, such as clus- 
ter age (Table HI), mid-IR luminosities (Table [3]), extinc- 
tion (Table EI), PAH strength (Table 2]), hardness of the 
radiation field (Table [5]) or star formation rate (Table [9]). 
However, no statistically convincing correlation between 
these quantities has been be found. 

4.6. Cluster Luminosities and Star Formation Rates 

The mid-IR continuum fluxes can be used to estimate 
the total infrared luminosity of starburst systems (in 
which LiR « Lboi)- In their Section 4.2 iBrandl et al.l 
((2OO6 ) discuss a method based on the 15/im and 30/im 
continuum fluxes. Adopting their empirical relation and 
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Fig. 7. — Comparison of the 14 — 37itm continuum slopes for 
three selected objects: peaks 1 and 5 with the most shallow and 
steep slopes in our sample of Ant ennae clusters , and the average 
starburst template spectrum from IBrandl et all ISiOOS) . The spec- 
tra are plotted in units of uF,, to enhance the differences, and have 
been normalized to unity at I Sjtm. 

a distance of 22 Mpc (jSchweizer et al.l l2008[ ) we get 
LiR = 2.067 • Fi5 + 5.324 • F30 in units of IO^Lq. We 
note that this empirical relation was derived for global 
starbursts, where the temperature components may dif- 
ferent from the one in and around more compact SSCs. 
However, the good agreement in Fig.[7]between the slopes 
of the various systems indicates that the estimate works 
reasonably well. 

In Table [H] we list the mid-IR derived cluster lumi- 
nosities. The infrared luminosities of the eight com- 
pact regions already add up to iiR — 3.82 • 10"'^° Lq. 
The total infrared luminosity of the Antennae, mea- 
sured from the f our IRAS bands at our distance, is 
LiR = 7.2-10^° LgllS anders et alll2003f ). Hence, the eight 
individual regions studied here account for slightly more 
t han half of t he tot al luminosity of the system. 

iKennicuttI (|1998D has shown that the 8 - lOOO/xm in- 
frared luminosity Lir of starbursts is a good measure of 
the SFR as given by 

SFR [Moyr^i] = 4.5 • IO-'^'^Lir [erg s'^]. 

We no te that the SFR conversion given by IKennicuttI 
(|1998[ ) strictly applies only to dusty starbursts in the con- 
tinuous star formation approximation, with ages of order 
10 — 100 Myr. Our clusters are signiflcantly younger than 
that, so in the continuous star formation approximation 
(which is a reasonable one, given the more or less con- 
tinuous range of cluster ages, as seen in Fig. |S|) the IR 
luminosity per unit mass of stars formed will be some- 
what lower. In other words, for regions younger than 
10 Myr the IKennicuttI (| l998l ) relation will systematically 
overestimate the SFRs. 

Our estimated SFRs for the individual clusters are 
listed in Table [H Co-adding the eight regions we get 
a total SFR w 6.6M0yr~^. Give n the compa c t size s 
of the regions our results agree with lWang eTal] (|2004'), 
who derived from Spitzer-ffiAC 3 — 8/xm observations 
that the rates of star formation in the active regions are 
as high as those seen in starburst and some ultraluminous 
infrared galaxies on a "per unit mass" basis. 

As a cross-check we compare our SFRs to an in- 
dependent estimate by dividing the stellar masses of 
the clusters from the literature by their ages. For 
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TABLE 9 

Derived Cluster Properties 





^IR 


or tx 


1(112] 


1(112] 


iVH^xl2 ) 


IVi yri2 ) 


Position 


[109 Lq] 


[Mo/yr] 


[K] 


[K] 


[IO^Mq] 


[106 M0] 


nuc 4038 


3.67 


0.63 


329 


327 


1.78 


2.24 


nuc 4039 


1.86 


0.33 


378 


334 


3.97 


3.56 


peak 1 


11.41 


1.97 




302 




3.09 


peak 2 


10.40 


1.81 


288 


267 


2.82 


3.20 


peak 3 


4.35 


0.74 


311 


300 


1.97 


3.39 


peak 4 


1.29 


0.22 




374 




0.22 


peak 5 


3.86 


0.66 


296 


309 


3.67 


4.24 


peak 6 


1.31 


0.22 




360 




0.30 



Note. — For uncertainties in the H2 mass and temperature estimates see Sec- 
tion fX^ The S(3)flux from peak 1 was well detected but is uncertain due to a data 
artifact. 

^ Luminosity estimate foUowing lBrandl et all | |2006|) ^ from the 'lores' S(3) and S(2) 
line fluxes'^ from the 'hires' S(2) and S(l) line fluxes'^ from the 'lores' S(3) line fluxes 
and temperatures'^ from the 'hires' S(2) line fluxes and temperatures 



peak 1 iMengel et al.l (12 001') determined 3 • 10^ Mq and 
iSniiders fc van der Wer f (2007) fomid 1.1 - 1.2- IO'^Mq. 
For peak 2 iMengel et al. ( 2001) found 1.6 • 10^ Mg 
and ISniiders &: van der W erf (2007) measured 1.2 — 1.7 • 
10^ Mq. For average cluster masses of 2 • 10^ M© and 
1.5 • 10^ Mq for peaks 1 and 2, and ages of 2.5 and 
3.0 Myr (Table [8]) we estimate (continuous) star forma- 
tion rates of 0.8 and 0.5MQyr~-'^, respectively. These 
rates are roughly a factor of three below the SFRs de- 
rived from the mid-IR fluxes. However, we emphasize the 
substantial systematic uncerta inties in these estimates. 
ISniiders k van der Werl (120071) assumed a Sal peter IMF 
from 0.1 - 100 Mq, while IMengel et al.l (|200lD assume a 
lower mass cutoff of IMq, corresponding to a 2.6 times 
smaller cluster mass. Using the same IMF would lead to 
mass differences between those studies of factors of 6.5 
and 2.9 for peaks 1 and 2, respectively. 

The sum of the SFRs around the two nuclei and the six 
infrared-brightest regions (4.8 ■ 10^ pc^ in total) as listed 
in table [5] is 6.6M0yr~^. As stated above, our method 
will likely provide an overestimate of the true SFR. On 
the other hand, star formation will also occur outside 
our regions at a reduced rate. Hence, we consider our 
estimate of 6.6MQyr~^, within 50%, as representative 
for the SFR in the Antennae galaxies^ 

We also note that , Zhang, F all & Whitmord (|200H ) de- 
termined a rather lar ge SFR of 26 Mq y r~^ for a distance 
of 22 Mpc. Using the iKennicuttI (|1998[ ) relation between 
SFR and Lm this would correspond to a more than two 
times higher total infrared lumino sity than the one de- 
rived from the IRAS bands (j Sanders et al.l [20CK^1 . and 
would move the Antennae to t he class of LIRGs. Hence, 
we con sider the SFR derived bv lZhang. Fall fc Whitmorg 
((2OOII) an overestimate. 

Recently, ISmith et al.l (|2007c') investigated Spitzer 
photometry of a sample of 35 tidally distorted pre-merger 
interacting galaxy pairs selected from the Arp atlas. 
These authors found a very modest enhancement of the 
SFRs at early stages of ~ IMq yr~^ on average. Al- 
though tlie_^Antennae is still in an early stage of inter- 
action (jMihos fc Hernquist 1996) it shows an enhanced 
SFR, well above the sum of typical star formation in two 
normal s piral galaxies. W ith its ample amount of molec- 
ular gas (|Gao et al.ll2001[ ) the system is likely to become 
an even more spectacular starburst after one or two or- 



bits. 

4.7. PAH Strength and Aperture Effects 

PAHs are considered the most efficient species for 
stochas tic, photoelectric heat ing by UV photons in 
PDRs ()Bakes fc TielensI I1994D . They are usually a 
good tracer of starburs t activity in a statistical sense 
(^e.g.. iBrandl et al1l2006l . and references therein). How- 
ever, numerou s stud ies over the pas t two decade s (e.g . 
iGeballe et af] fl989t iCesarskv et all Il996l: ItVm] [1993) 
have shown that intense UV fields can lead to the gradual 
destruction of PAH molecules. A clear anti-correlation 
between the 11.3/im PAH strength and the radiation field 
across the central regio n of NGC 5253 was recently ob- 
served bv lBeirao et al.1 12006). We investigate this effect 
for the Antennae in two different ways: via a pixel-by- 
pixel analysis within the spectral map and for the eight 
'hires' spectra. 
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Fig. 8. — The relative strengtli of the 11.3;im PAH feature (a 
normahzed "line"-to-continuum ratio) versus the hardness of the 
radiation field as expressed by the line ratio of [S IV] / [No II] . Each 
data point corresponds to one pixel in the spectral map. Only those 
42 pixels with a signal of greater than four standard deviations in 
the S IV map have been included. 

The former is shown in Fig. |5| where the 11.3/xm 
PAH-to-continuum ratio is plotted versus the hardness 
of the radiation field as expressed by the line ratio of 
[S IV] / [Ne II]. Each data point in Fig. [8] corresponds to 
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a spatial pixel in our spectral map for which the signal-to- 
noise at [S IV] is greater than four standard deviations. 
(See section [531 concerning the contamination of [Ne II] 
by the 12.7/im PAH emission in 'lores' mode). Fig. [S] 
reveals a trend that regions with harder radiation fields 
show relatively weaker PAH emission. 

We note that the strength of the [S IV] relative to 
the [Ne II] line does not only depend on the hardness 
of the radiation field, but is sensitive to the gas density 
as well. Given the critical densities of 3.7 • 10* c m~^ and 
6.1-1 0^ cm-3 for [S IV] and [Ne II], respectively ()Tielensl 
l2005l and references therein), it is expected that den- 
sities above lO'^ cm~^ will decrease the [S IV] / [Ne II] 
ratio. Model calculations with CLOUDY show that an 
increase in density from 10* to 10^ cm~'^ can cause a de- 
crease in [S IV] / [Ne II] by a factor of three, contributing 
significantly to the scatter in Fig. [51 

To quantify the strength of the radiation field we use 
the lin e fluxes of [Ne II] and [Ne III] and the parametriza- 
tion of iBeirao etall ()2006l ): 



(-F'lNc II] 12.8^111 + ^'[NcIIIJlS.eAim) X 



^'[NeIII]15.6/ir 
-P'[NcII]12.8Am: 



The first factor, Ft 



[NcII]12.8/im 



[NoIII]15.6/im: 



IS a mea- 
sure of the intensity of the field, assuming that all Neon 
exists ii^one of the two ionization states. The second 
factor, jg a measure of the radiation field 

hardness. We define the product of intensity and hard- 
ness as the "strength" of the radiation field. 
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Fig. 9. — Tlie relative strength of all PAH features within the 
10 — 20fj,m range versus the strength of the radiation field. 



In Fig.iniwe plot the strength of the PAH features, mea- 
sured from the 'hires' spectra, against the strength of the 
radiation field. On the ordinate we plot the average PAH 
strength calculated by co-adding all the PAH fluxes listed 
in Tableland normalizing them to the 15/im continuum 
flux (Table [3]) . There is a clear linear anti-correlation 
between the strength of the PAH emission and the ra- 
diation field. The only outlier is the inactive nucleus of 
NGC4038. 

Both rather independent methods and data sets pro- 
vide statistical evidence that the strength and hardness 
of the radiation field affect the PAH emission. However, 
we cannot be certain that this effect is actually due to a 
destruction of PAH molecules in stronger radiation fields. 



The sources with the strongest PAH "equivalent width" 
are peaks 4 and 6. Both are located in less dense re- 
gions. These clusters are also slightly older and have 
had more time to shape the surrounding ISM via stel- 
lar winds and SNe. If the size of the emitting region is 
radiation-bounded these regions may just have "grown" 
more complex Hii region/PDR interfaces from where the 
PAH emission originates. In other words, the PDR sur- 
face area is higher relative to the thermal emission from 
dust in the denser regions, which produces the underly- 
ing continuum. 

Unfortunately, with the wide ffiS-SH slit width, cor- 
responding to 500 pc, we cannot resolve the Hii re- 
gion/PDR interfaces. However, some of the clusters were 
recently observed by [Sniiders et"all (|2006[ ) with VISIR 
on ESO's VLT using a sub-arcsecond slit correspond- 
ing to only 75 pc at our distance of 22 Mpc. Compar- 
ing our IHS-SH spectra of peaks 1 a nd 2 to their VLT- 
VISIR spectra, [Sniiders et all (|2006[ ) found that, while 
the continuum fluxes observed with both instruments are 
quite similar, the wider IRS slit detects much stronger 
PAH emission. The individual clusters are by far the 
most luminous sources on projected scales of a few hun- 
dred parsecs, and we may observe the profound effect 
of the OB clu s ters o n their surrounding ISM. However, 
iSni iders et al.[ ([2006[ ) concluded that a large fraction of 
the PAH emission cannot be directly associated with the 
SSCs but originates from an extended region around the 
cluster. In fact, if the space density of young clusters 
near peaks 1 and 2 is significantly higher than elsewhere 
it could be that the extended PAH emission, which is only 
picked up by the IRS, is unrelated to the main clusters 
and externally excited by smaller clusters in the imme- 
diate vicinity of peaks 1 and 2. At this point we cannot 
distinguish between a giant Hii region of a hundred or 
more parsecs in diameter and a locally enhanced density 
of OB clusters. A significant improvement will have to 
wait for the next generation of extremely large telescopes 
(ELTs). 

4.8. (In-) Variability of the PAH Spectrum 

We have also investigated possible variations of the 
PAH spectrum. While - to first order - the strengths 
of the PA H features are obse rved to scale with each 
other fe.g.. [Smith et al.|[2007bl Fig. 6), the PAH spec- 
trum is expected to depend on the grain size distribu- 
tion (Draine fc Lf 2001, 2007): smaller PAHs are more 
likely to emit at shorter wavelengths while larger PAHs 
are more likely to produce stronger features at longer 
wavelengths. We utilize the longest baseline provided 
by the 'hires' spectra between strong PAH features and 
selected the 11.3/xm feature and the 17fim PAH com- 
plex^'^. We find an average 17/im / 11.3/im PAH ratio of 
0.4 (Table[3]) with some scatter but no significant correla- 
tion with any othe r observed quantity. From the SINGS 
sample of galaxies ISmith et al.[ ()2007bl ) have found that 
the 17/im / 11.3/im PAH ratio increases with increasing 
metallicity. Assumin g solar abundance for the Antennae 
(jBastian et al.|[2b06[ ) our PAH ratio of 0.4 is consistent 

The 17Atm P AH complex includ es the 16.4^ m feature 
IllVIoutou et al.ll2000l) . the 17.4/im feature llSturm et al. 2000), and 
a very broad feature at I7 .lfim, which cont ains over 80% of the 
PAH flux of that complex llSmith et al.ll200l '). 



16 



Brandl et al. 



with the measurements of ISmith et alj ()2007b . Fig. 16) 
although somewhat on the low side. 

Another cause of potential variations of the PAH 
spectrum r aay be ionization effects. Several au- 
thors (e.g., IVe rstractc ct al.J Il996t IVermeii et all l2002t 
iForster Schreiber et al.. .2003? ) have reported on varia- 
tions in the PAH spectrum, namely that the C-C stretch- 
ing modes at 6.2/.tm and 7.7/.tm are stronger in ionized 
PAHs, relative to the bending mode at 11.3/im - by fac- 
tors of up to two. Here we use the ratio of the 8.6/im C-H 
in-plane bending mode to the 11.3/im C-H out-of-plane 
bending mode, which is a.lso sensitive to the charg e state 
of the PAH mol ecule (jHudgins fc Allamandolai Il995t 
iJoblin et al.l |1996[ ). We have chosen these two for a 
direct comparison because of their proximity in wave- 
length (same beam size), similar susceptibility to dust 
extinction, and becau se both are o bserved through the 
same IRS (sub-)slit. ISmith et all (|2007bD investigated 
the 7.7/im / 11.3/im PAH ratio as a function of the hard- 
ness of the radiation field and found no significant cor- 
relation for the Hii-nuclei in the SINGS galaxies sample 
(their Fig. 14). 

Fig. [4] already indicates that, to first order, both PAH 
maps trace each other quite well. Thus we have com- 
puted the ratio on a pixel-by-pixel basis between the 
two 8.6/im / 11.3/im PAH maps. The map is shown in 
Fig.[Tn]at two times lower resolution to reduce the noise. 
Regions with relatively stronger 8.6/im emission appear 
brighter, regions with relatively stronger 11.3/im emis- 
sion appear darker in this ratio map. While the emission 
from the nucleus of NGC 4038 appears in average "gray" 
the peaks 3,5, and 6 appear darker and have values about 
20% above the mean. On the other hand, peaks 1 and 2, 
the most massive clusters, appear brighter by about the 
same amount. These variations are not correlated with 
extinction effects, which would affect the values in the op- 
posite direction. However, we emphasize that the noise, 
in particular in the 8.6/im map (Fig. [5]), is quite signif- 
icant. Although these variations may provide weak evi- 
dence for changes in the PAH spectrum on cluster scales, 
a conclusive analysis requires better resolution and lower 
noise. 

4.9. Strength of the H2 Emission 

Molecular hydrogen is the most abundant molecule 
in the Universe. It plays a central role in star forma- 
tion, not only as the major ingredient to build up a 
star but also as the coolant to permit an isothermal 
collapse of the gas cloud. The density and tempera- 
ture of the H2 molecule are of utmost importance for 
the processes in starbursts. However, the H2 molecule 
is symmetric, has no dipole moment, and the mid-IR 
lines originate from quadrupolar rotational transitions. 
Consequently, the lines are intrinsically weak and usu- 
ally hard to dete ct. Since the fi rst study of the S(2) 
line in Orion by B eck et all (Il979f) and the detecti on of 
the S(0) fine in NGC 6946 bv lValentiin eTall (fl99l . the 
mid-IR rotational lines in starbursts have been studied 
by numerous groups. Most relevant in our con text are 
the ISO-SWS studies bv'Rig opoulou et al.l (|2002'), the re- 
view by Habart et al. (2005), and t he Spitzer-JfiS studies 
of ULIRGs by Higdon c t alJ pOOl . and SINGS galaxies 
by Roussel et al.. (2007, ). 

The H2 lines in the Antennae are amongst the strongest 




Fig. 10. — IRS spectral map of the 8.6/im / 11.3/im PAH ratio. 
The median ratio has been normalized to unity, which is colored 
in mean gray, and the range extends over a factor of two to both 
sides. Brighter regions represent relatively stronger 8.6/tm emis- 
sion, darker regions weaker 8.6/im emission. Field size and ori- 
entation are as in Fig. U] but the resolution has been reduced by 
rebinning by a factor of two to suppress the noise. 

detected lines (Fig. We have calculated the total, 
spatially integrated H2 S(2) and S(3) fine fluxes by sum- 
ming up all pixels where the signal-to-noise was greater 
than three standard deviations in the spectral maps and 
compared them to the sum of the regions listed in Ta- 
bled In the S(3) line essentially all of the measured flux 
of 6.1 • 10~^° W cm"^ originates from our designated 
regions. The emission in the S(2) line is slightly more 
extended and the spatially integrated measure picks up 
27% more flux than listed in Tabled We conclude that 
most of the total amount of warm H2 over the central An- 
tennae system is concentrated around our eight sources. 

Fig. E] shows a false-color image of the IRS maps 
with the H2 S(3) and S(2) lines in blue/green and red. 
Both nuclei and the emission from the active overlap re- 
gion are clearly visible. The map shows that the largest 
portion (approximately 45%) of the total H2 S(3) emis- 
sion comes from a region around the southern nucleus. 
The strongest H2 emitter in the active overlap region is 
peak 5. 

We note that our integrated H2 S(3) line flux of 
6.1 • 10^^° Wcm~^ is about five t imes below t he S(3 ) line 
flux of 3.3 ■ 10-19 Wcm-2 which iHaas eTal] (|2005[) have 
found. In Fig. [TT] we overlaid the contours of their ISO 
H2 S(3) intensity map for comparison. While our peak 5 
is close to their southern peak of the contour map, the 
northern peak of the contour map has no luminous coun- 
terpart in our map. Our closest peak is peak 6, which 
contains more than an order of magnitude less H2 mass 
than peak 5. We have also checked the IRS-LL map, 
which includes the H2S(1) emission line, albeit at two 
times lower spatial resolution, and confirmed the strong 
drop in H2 emission northward o f peak 3 . Altogether, we 
cannot confirm the claim by Haa s et all (|2005) of strong 
H2 emission from a region extending farther to the north 
of the overlap region, well beyond the most active area. 
These differences in distribution and total line fluxes have 
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Fig. 11. — False-color IRS spectral map of the H2 S(3) line in 
the green and blue channel, and the S(2) line in the red channel in 
square-root scaling. Bluer colors mean warmer H2. North is up. 
The two nuclei and the emission from the clusters in the overlap 
region are clearly visible. The yellow contour lines i ndicate the 
ISOC AM-CVF measurements of the H2 S(3) line by IHaas et al.l 
lHOOl); the levels are at 39%, 49%, 58%, 68%. 78%, 84% 90%, a nd 
96% of the maximum intensity (see Fig. 2 in IHaas et al.l (120051 )). 

significant impact on the derived H2 masses for the An- 
tennae system. 

At the distance of 22 Mpc, our S(3) line flux corre- 
sponds to a line luminosi ty of 9.2 • 10^ Lp, . Usi ng the far- 
IR luminosity derived bv lSanders et al.1 (|2003f) we calcu- 
late the ratio to be L(H2)/-Lfir ~ 1-6 • 10~^. This value 
is only about two times higher than the one for M82 and 
not atypical for luminous starbursts and ULIRGs (cf. 
IHaas et al.l[2005l Fig. 3). 

Here we also investigate with what other observ- 
ables the strong H2 e mission might be correlated. 
iRigopoulou et alT ()2002[ ) found a good correlation be- 
tween the 7.7/Ltm PAH emission and the H2 S(l) line lu- 
minosity in ULIRGs. A tight correl ation between H2 and 
PAH emission was also reported bv lRoussel et al.l ()2007f ) 
for the SINGS galaxies (their Figs. 8c and 9). In fact, 
they found that "among the measured dust and gas ob- 
servables, PAH emission provides the tightest correlation 
with H2". This sounds plausible since both PAHs and 
H2 reside preferentially at the outer skin of PDRs, and 
are both predominantly exited by FUV photons. How- 
ever, while PAH emission is almost exclusively excited by 
UV photons, H2 can be excited by shocks or collisionally 
excited if the critical density of typically ~ lO'^ cm~'^ 
is exceeded (Section I4.12p . This may well be given lo- 
cally for the embed ded SSCs where e lectron densities 
up to 4 ■ 10"^ cm"'^ (I Schulz et all l2007l ) or even higher 
(jMirabel et al.|[T998[) have been measured. 

In Fig. [12] we investigate a possible correlation between 
the H2 line strength and the PAH strength. The com- 
parison has been done for the eight sources observed in 
'hires' mode and for the spectral maps on a pixel-by-pixel 
basis. For both cases we do not find a tight correlation. 
In fact, the apparently uncorrelated scatter in Fig. [T2lis 
remarkable, although not completely unexpected given 
the different morphology between the PAH and H2 maps 
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Fig. 12. — The strength of the H2 lines plotted versus the 
strength of the PAH emission features. The figure contains two 
sets of data: the large diamonds correspond to the eight 'hires' 
positions for which all PAH features longward of lO^tm have been 
added; the dots refer to pixels within our 'lores' spectral maps for 
which the signal is at least 2.5o-. In the case of the 'lores' data the 
PAH strength is the sum of the 8.6^m and 11.3/jm maps and the 
H2 strength has been derived from the S(2) and S(3) lines. The 
pixel fluxes have been scaled to match the 'hires' data points. 

in Fig. H We note that iRoussel et"al] (|2007f) base their 
PAH index on the 7.7/im feature while we use different 
PAH features. However, since we do not observe any sig- 
nificant variations within the PAH spectrum (section l475|) 
we assume that both measures are equivalent. 

There are two possible explanations for the difference 
between the previously observed PAH-H2 correlation and 
our results. First, the sources of strong H2 emission in 
the Antennae are compact regions of a few hundred par- 
sees at most. Larger scatter aroun d the average galacti c 
value is therefore not unexpected (jRoussel et al.l 120071 ). 
Second, we cannot rule out contributions from local 
shocks caused by supernovae. The latter scenario will 
be discussed in Section B.121 

4.10. H2 Temperatures 

From a sample of 77 ultra-lumin ous infrared g alaxies 
(ULIRGs) observed with the IRS iHigdon et alT (|2006h 
measured a mean temperature of the warm H2 gas of 
336 K. For the more quiesc ent galaxies in the SINGS 
sample, IRoussel et al.l (|2007l ) determined a median tem- 
perature of only 154 K. However, the slit apertures of 
the SINGS observations sample the circum-nuclear re- 
gions; often these spectra have contributions from mul- 
tiple emitting sources, and presumably regions with a 
range of dens ities and t emperatures. For the Anten- 
nae galaxies K unze et al. (1996) deri ved a te mperature of 
405 K from ISO-SWS data, whereas [Fischer et al.l (|l996i ) 
used ISO-LWS (while using different tracers) and derived 
a temperature for the PDR of 200 ± 60 K. We note that 
the common decomposition in warm and cold H2 compo- 
nents may not be unique and the absolute temperatures 
not necessarily physically distinct. However, they allow 
useful qualitative comparisons between different systems. 

Our temperature estimates for the different regions are 
listed in Table [HI We emphasize that the two indepen- 
dent temperature estimates for the IR peaks which were 
derived from two different sets of lines agree remarkably 
well, to within a few percent. This is likely because the 
individual spectra are dominated by single sources, which 
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almost completely thermalize the warm H2 . Such a good 
agreement between temperature estimates from the S(l), 
S(2) and S(3) Hnes woul d not be expected if th e H2 were 
primarily shock-excited (jAppleton et al.ir2006[ ). 

To compare our estimates with the values from the 
literature we "simulated" a larger slit aperture (e.g., of 
ISO-SWS) by calculating a luminosity-weighted mean 
temperature from the S(2) and S(l) lines for the clus- 
ters peak 1, 2, 3, 5, and 6 within the overlap region. 
This yields a mean warm H2 temperature of 302 ± 26K, 
which lies sig i iifican tly below the temperature estimate of 
iKunze et al.l (|1996f l for the Antennae, but only 10% be- 
low t he average temperature of ULlRGs (jHigdon et al.l 
|2006[ ). The quoted uncertainty on the temperature is 
the 1 — a standard deviation between the temperature 
estimates in Table [9l and includes the temperature vari- 
ations between clusters as well as the independent mea- 
surements from different IRS modules. Hence, we con- 
sider this uncertainty as a good estimate of the total 
systematic uncertainties. 

We also investigated the possible dependency of the 
temperature of the warm H2 component on other ob- 
servables. Fig. [T3l shows the relative strength of the PAH 
emission, the hardness of the radiation field, and the star 
formation rate as a function of the H2 temperature. All 
three plots show interesting correlations. 

For each of the eight regions we have computed the 
sum of all PAH features between 10 — 20/im and normal- 
ized them to the continuum flux at 15/im. The relative 
PAH strength correlates well with the H2 temperature. 
It appears likely that those PDRs, which are stronger 
PAH emitters relative to the continuum emission, have 
warmer molecular gas. On the other hand, the hardness 
of the radiation field (center plot in Fig. [T3|) , expressed 
in terms of [Ne III] / [Ne II] , anti-correlates with the H2 
temperature: regions of higher temperature show a softer 
radiation field. The general decrease toward hotter H2 is 
also observed individually for both line fiuxes of [Ne II] 
and [Ne III], but is steeper for the latter, resulting in the 
decrease of radiation hardness in hotter regions. 

These findings may, at first, be surprising as one might 
have guessed that a harder radiation field is more ener- 
getic and will thus produce a hotter environment. Our fa- 
vored explanation is that the older, more evolved clusters 
had more time to shape the surrounding ISM via stel- 
lar winds and supernovae, and the surrounding, clumpy 
dust and gas clouds have evolved into larger entities. The 
UV photons from the cluster can penetrate farther into 
the cloud and provide a more uniform temperature dis- 
tribution over a larger area. Very young clusters may 
have hotter gas and dust near their inner shell but are 
colder at their well shielded outskirts. This is mainly a 
geometrical argument, similar to the scenario discussed 
in Section 14.71 We also note that the top plot shows 
the relative PAH strength, normalized by the continuum 
fiux at 15/im. The absolute PAH fluxes, however, show a 
similar trend as the Neon lines, namely stronger emission 
from regions of colder gas. This indicates that the hotter 
regions are more eflScient PAH emitters. 

Finally, the bottom plot in Fig. [13] shows the star for- 
mation rate (Table [9|) versus the H2 temperature. Al- 
though less significant than in the upper two plots there 
is a clear trend of lower SFR for hotter regions. This 
is consistent with the previous anti-correlation between 
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Fig. 13. — Dependency of three important parameters on the 
temperature of the H2 gas, calculated from the S(2) and S(l) lines. 
Top; The relative strength of the PAH emission (given by the flux 
sum of all PAH features between 10 — 20/jm normalized to the 
continuum flux at IS/^m). Center: The hardness of the radiation 
field (as given by the line flux ratio of [Ne III] / [Ne II]). One outlier, 
the nucleus of NGC 4038, which shows a very low [Ne III] / [Ne II] 
ratio, is not included. Bottom: The star formation rate as listed in 
Tabled The dotted lines are linear fits but only intended to guide 
the eye. 

both hardness and line fiuxes - which are expected to 
be stronger in regions of recent star formation, and the 
H2 temperature. Since the SFRs have been derived inde- 
pendently the good agreement with our line diagnostics 
underlines the reliability of the results. 

4.11. H2 Masses 

The fraction of the warm to total molecular hydrogen 
mass is of particular interest as it provides information 
on the dominant excitation mechanism and can be com- 
pared to other galaxies. The total molecular hydrogen 
mass is dominated by the cold component, which is de- 
rived from ^^CO measurements and then converted to 
H2 using the standard Galactic conver sion factor. 

For starburst-dominated galaxies iRigopoulou et all 
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(|2002f) found that the warm H2 (T ~ 150 K) accounts 
for about 1 — 10% of the total molecular hydrogen mass. 
In Seyfert galaxie s this fraction can be higher (2 — 35% 
iRigopoulou et all 120021. Recently, ISchulz et all (|2007f ) 
reported that the gas properties in the Antennae do not 
deviate significantly from the f^n-^lico ratio found for 
the disk of the Milky Way, an d the convers i on fro m CO 
is thus reasonably accurate. IWilson et al.l (|2000[ ) used 
a CO to H2 conversion factor of 3 • 10^° H2Cm~^, and 
derived a total H2 mass of 7.1 • 10^ Mq, ado pted to our 
distance of 22 Mpc. Earlier estimates by St anford et al.l 
(|1990t ) yielded a total H2 mass of 1.3 • 10^ M©. 

The derived masses of warm H2 are listed in Table [9l 
We note that the two independent estimates from dif- 
ferent ffiS modules, H2 lines and observing modes agree 
within 15%. Adding up the warm molecular gas masses 
for the 'hires' positions in Table [9] yields ][,f^9^°"'' — 
2.0 • 10^ M0. If we also account for the more extended 
"missed" gas (approximately 25% in the S(2) line - Sec- 
tion[49l) we get M|g °' = 2.5-IO'^Mq. Ba sed on the total 
mass estimate from IWilson et al.l (|2000( ) we measure a 
fraction of warm-to-total molecular gas mass in the An- 
tennae of 0.35%. This ratio needs to be compared to 
the ex t reme case of NGC 6 240 where Ivan der Werf et aD 
(|1993l ): lArmus et al.l (|2006f ) found about 15% of the total 
molecular gas in its warm phase. However, the Antennae 
appears to have, for its luminosity, an average fraction 
of warm-to-total molecular gas mass. 

These estimates can also be made for individual 
regions. The supergiant molecular clouds labelled 
S GMC 1, SGMC 2 and (SGMC 3 + SGMC 4 + SGMC 5) 
bv IWilson et al.l (|2000( l correspond approximately to our 
peak 3, peak 5, and (peak 1 + peak 2), respectively. 
Moving both systems to the same distance, we calculate 
fractions of warm to total gas masses of 0.4%, 0.8%, and 
0.6% for the three regions, respectively. The more local- 
ized ratios are somewhat higher than the galactic aver- 
age, which is not unexpected since the cold H2 is likely to 
cover a larger volume. Similar results are obtained from a 
comparison with the regions 'south c lump' , 'west clump', 
and 'east clump' in IStanford et al] (|1990D . which corre- 
spond approximately to our (peak 1 -|- peak 2), peak 5, 
and peak 3, respectively. 

4.12. H2 Excitation Mechanisms 

All the observable lines at mid-IR wavelengths origi- 
nate from quadrupolar rotational transitions. The ro- 
tational levels of H2 can be populated in va rious ways: 
H2 molecules may form in an excited state ()Takahashil 
12001:): X-ray photons from energetic sources like AGN 
or supernova remnants can ionize and heat the gas, 
which excites the H2 molecules via collisions with elec- 
trons and /or hydroge n atoms ()Lepp fc McCravl 119831 : 
iDraine fc Wood !j'l992^: FUV p hotons can pump the H2 
molecules in ex cited states (Bl ack fc van D ishocck 1987; 
iHoUenbach fc T iclcns 1997); a nd excitation via shocks i n 
the ISM (iShtll fc Hollenbachl[T97at IPraine et al.lll983h . 
In starburst environments we can neglect the former two 
scenarios and concentrate on UV pumping and shock ex- 
citation. 

UV pumping requires photons with energies E^y = 
6 — 13.6 eV to excite the H2 molecules, which cas- 
cade via vibration-rotational transitions, producing char- 



acteristic fluorescent spectra. Shocks are expected 
to occur on various scales. On small scales, shocks 
from stel lar outflows and superno va remnants can ex- 
cite H2 (jShull fc HollenbachI [19781) . On larger, galac- 
tic scales streaming motions and cloud-cloud collisions, 
in particular in merging systems have been proposed 
(Ivan der Werf e t al. 199 3[lGilbert et al.ll2000l:lllaas et all 
l2005HArmu~ al.„2006i : lAppleton et al.ll2006D ~ 

Both FUV pumping and shock excitation have been 
shown to play a dominant rol e in the excitation of H2 in 
starburst enyironni ents (e.g.. iMoorwood fc Olival [19941 : 
IRoussel et al.ll2007l and references therein). Pure fluo- 
rescent spectra can usually be distinguished from ther- 
malized, shock heated spectra from the relative line 
strengths. However, depending on the critical densi- 
ties, collisional de-excitation can also become impor- 
tant for UV pumped H2 states, leading to thermalized 
lower levels. Most of the emission of warm H2 is likely 
to come from dense PDRs with densities > lO'^cm"^, 
and we expect that the observed S(0) through S(3) 
lines are thermalized (jBurton et al.l il992'). The near- 
IR rotation-vibration transitions 2-1 S(l) at 2.248/im 
and 1-0 S(l) at 2.122/im can be used t o discriminate 
betw een shocks and UV excitation (e.g., iTakami et al.l 
|2000() . Their ratio is typically 0.1 in shocked regions 
(with the gas in LTE at T ^ 2 ■ lO'^ K), but ^ 0.6 fo r 
"pure" UV fluorescence (|Black fc van DishoecH I1987D . 
However, in very dense PDRs (~ 10^ cm~'^) even the 
near-IR rotation-vibrational lines are not a good dis- 
criminant anymore as thermal collisions can transfer the 
lower-level (z/ < 2) populations toward that in LTE, and 
will resemble the ratio obser ved in shocked regions (e.g., 
iSternberg fc Dalgarnolll989[ ). 

The observa tional evidence in the Antenna e so far 
is confusing. 'Sni iders fc van der Werj (|2007f) derived 
2-1 S(l) / l-0S(lj values between 0.15 and 0.34 for 
peaks 1, 2, and 5. These values are in the "grey 
zone" between the model values for UV pumping and 
shocks. However, thermal collisions may lower these 
ratios, and the presence of high iz-level lines {ly = 
4, 5, 6 , 7) suggests that fluorescence pla ys an important 
role (ISniiders fc van der Werj 120071 ). IWhitmore et all 



()2005[) derived a velocity dispersion between star clus- 
ters of < lOkms^^ from STIS spectroscopy and con- 
cluded that high-velocity cloud-cloud collisions cannot 
be the dominati ng starburst mode i n the Antennae. On 
the other hand, IWilson et al.l (|2000D found evidence for 
clou d-cloud collis i ons n ear the strongest mid-IR peak, 
but ISchulz et al.l ()2007f ) concluded that massive shock 
heating, either due to supernova remnants or due to 
cloud collisions, is not likely the dominant heating source 
for the SGMCs in the overlap region. 

A starburst system of particular interest for compar- 
ison may be Stephan's Quintet (NGC 7317/7318b)- a 
compact group of four strong ly interacting ga l axies with 
a likely foreground galaxy. lAppleton et al.l ()2006D re- 
ported on the discovery of an almost "pure H2 line spec- 
trum" with unusually strong and broad H2 lines, resem- 
bling the spectra of shocked gas seen in Galactic super- 
nova remnants. The total mass of warm H2 was deter- 
mined to be 3.4 • 10^ Mq, about 40% more tha n what 
we measured in the Antennae. lAppleton et al.l ()2006D 
attributed the powerful H2 emission from Stephan's 
Quintet to a large-scale shock wave, resulting from a 
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high-velocity galaxy collision. The evidence for this 
claim is manifold: the unusually broad line width of 
~ 870kms"\ the relative strength of the H2 S(0) - S(5) 
lines as indicated in the excitation diagram (their Fig. 3), 
and the absence of PAH-dust features and very low ex- 
citation ionized gas tracers. 

So far, the strongest support for large-scale shock exci- 
tatio n in the Antennae came from ISO data (iHaas et al.l 
|2005[) . which showed exceptionally strong H2 S(3) emis- 
sion from the overlap region. In addition, the detected 
emission appears also spatia lly displaced from the known 
starburst regions (Fig. [11]). iHaas et all (|2005| ) interpret 
these two findings as combined evidence for pre-starburst 
shocks that arise from neutral H i cloud-cloud collisions. 

As discussed in Section l4?9l ou r IRS obseryatiori s dis- 
agree with the measurements by IHaas et al.l ()2005D . We 
find that the H2 emission appears to be reasonably well 
correlated with the starburst activity in the southern 
part of the overlap region. The only exception may be the 
nucleus of NGC 4039 for which - to our knowledge - no 
published near-IR H2 flux me asure ments exist, but the 
spectra of iGilbert et~an (|2000D and lMengel et all (|200lD 
suggest that the ratio is very small and lies within the 
'shocked' regi me. Based on the n ear-IR lines of H2 1-0 
S(l) and Br^, iFischer et al.l (|1996f ) argued that the sur- 
face brightness at the nucleus of NGC 4039 is too high to 
be explained by UV excitation, and favoured C-shocks. 

In summary, we find no evidence for large-scale pre- 
starburst shocks due to Hi cloud-cloud collisions. We 
also note that we have seen none of the strong evi- 
dence for large-scale shocks, that has been reported for 
Stephan's Quintet bv lAppleton et all ()2006[ ). in the An- 
tennae. In the overlap region the H2 emission agrees well 
with the starburst activity, and leaves both FUV photons 
and local shocks, e.g. from supernovae, as origin for H2 
excitation. In the southern nucleus, from which 45% of 
the H2 S(3) is detected, shocks appear to be the domi- 
nant trigger. Further observations at higher resolution, 
e.g. by JWST, are necessary to get a clearer picture of 
the complex processes in this interacting system. 

5. SUMMARY 

We observed the Antennae galaxies (NGC 4038/39) 
with the Infrared Spectrograph on board of the Spitzer 
Space Telescope to study the properties of the ISM 
and the most luminous super-star clusters in this pro- 
totype merger system. We obtained low-resolution (i? ~ 
100) spectral maps of the entire central area and high- 
resolution {R ~ 600) spectra of six infrared-luminous 
regions and the two galactic nuclei. The high signal-to- 
noise spectra allow for a detailed study of fine-structure 
lines, PAHs, silicates and molecular hydrogen in the 
range between 5 — 38/im. 

Both nuclei of NGC 4038 and NGC 4039 are sur- 
rounded by regions of moderately active star formation. 
While the nucleus of NGC 4039 appears more active in 
terms of the [Ne III], [S IV] and H2S(3) emission lines, 
we can rule out the presence of an AGN in both nu- 
clei, based on the lack of [Ne V] emission and the low 
[O IV] /[Ne II] ratio. Generally, the older super star clus- 
ters are associated with softer radiation, as measured by 
the [Ne III]/[Ne II] ratio. The by far hardest and most 
luminous radiation originates from two compact clusters 
(our peaks 1 and2) in the southern part of the overlap re- 



gion. The tracers of softer radiation and PAH emission 
- typical for photo-dissociation regions ~ are spatially 
extended throughout and beyond the overlap region. 

From the slope of the continuum emission we found 
that peaks 1 and 2 have the highest dust temperatures 
while peaks 3 and 5 have the lowest. However, the dif- 
ferences are subtle and in agreement with the average 
spectral slope and diversity of starburst galaxies. The 
total infrared luminosity of our six IR peaks plus the 
two nuclei is Lir = 3.82 • lO"'^" Lq. From the IR luminosi- 
ties we derived the star formation rates for the individual 
regions, which range between 0.2 and 2 MQyr^^ with a 
total of 6.6 Mgyr-^ 

We found weak evidence for suppressed PAH emis- 
sion from regions with harder ([Ne III]/[Ne II]) emission. 
However, there is a clear trend that regions with stronger 
(i.e., harder and intenser) radiation fields show reduced 
PAH strength. Since the Hii regions and PDR interfaces 
cannot be spatially resolved at the distance of the Anten- 
nae, aperture effects play an important role in the inter- 
pretation of the results. We found no evidence for PAH 
grain size variations from the comparison of the 17/im 
and 11.3/zm PAH features. To first order, the PAH fea- 
tures trace each other very well. However, we find some 
evidence for spatial variations in the 8.6/im / 11.3/im ra- 
tio, which could be attributed to PAH ionization effects. 

The H2 lines are among the strongest emission lines 
in the Antennae. The emission is rather confined to the 
region around the nucleus of NGC 4039 and the south- 
ern, more active, part of the overlap region. We find a 
total line flux of 6.1 • lO^^o Wcm^^ in the S(3) li ne. This 
is a fa ctor five less than previous claims by iHaas et al.l 
(Unni) based on ISOCAM-CVF data. With an S(3) line 
luminosity of 9.2 • 10^ Lq the luminosity ratio between 
the warm molecular gas (as measured by the H2 lines) 
to the total FIR emission is 1.6 • 10~^, a value similar 
to that found in many starburst and ultra-luminous in- 
frared galaxies. We did not find the tight correlation 
be tween the H2 line and the PAH strength that was seen 
by iRoussel et al.l (|2007t l on larger scales in the sample 
of SINGS galaxies. Apart from the general "smoothing 
effect" when averaging over larger regions local shocks 
from supernovae are likely to contribute to the discrep- 
ancy. 

We derived a temperature of the warm H2, averaged 
over the individual star forming regions, of 302 ± 26 K. 
PDRs which have stronger PAH emission relative to the 
continuum emission have warmer molecular gas. Inter- 
estingly, the hardness of the radiation field anti-correlates 
with the H2 temperature. We explain this by the PDR 
geometry: older clusters had more time to shape their 
surroundings; the UV photons can penetrate farther into 
the cloud and provide a more uniform temperature dis- 
tribution, while young clusters may have hotter gas and 
dust near their inner shell but appear cooler at their well 
shielded outskirts. We calculated a total mass of warm 
H2 in the Antennae of 2.5 • 10^ Mq. The fraction of warm 
to total H2 gas mass is 0.35% globally and varies between 
0.4% and 0.8% for the individual supergiant molecular 
clouds. 

We find no evidence for large-scale pre-starburst shocks 
due to H I cloud-cloud collisions. Whether the H2 emis- 
sion is due to local shocks or UV pumping cannot be 
unambiguously determined from our data. In the over- 
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lap region the H2 emission agrees weh with the starburst 
activity, and leaves both FUV photons and local local 
shocks, e.g. from supernovae, as origin for H2 excitation. 
In the southern nucleus, from which 45% of the H2 S(3) 
is detected, shocks appear to be the dominant excitation 
mechanism. 
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TABLE 10 

H2 MOLECULAR CONSTANTS 



A [/im] 




9s 


J up 


9J 


Eup/k [K] 




- [s-i] 


28.219 


S(0) 


1 


2 


5 


509.9 


2.94 


■10-" 


17.035 


S(l) 


3 


3 


21 


1015.1 


4.76 


■ 10-10 


12.279 


S(2) 


1 


4 


9 


1681.7 


2.76 


■10-9 


9.665 


S(3) 


3 


5 


33 


2503.8 


9.84 


■10-9 



^ The upp er-level en ergies and Einstein A coefficients were 
taken from lRosenthal et al. ( 2000. ) and references therein. Qs 
is the statistical weight factor and gj is the level degeneracy 
of rotational level J. 

APPENDIX 

CALCULATION OF THE H2 TEMPERATURES AND MASSES 



The H2 molecule exists in two states: ortho-H2 with parallel nuclear spins (odd J), and para-H2 wit h anti-parallel 
spins (even J) . In local thermodynamic equilibrium (LTE) the ortho-to-para ratio is 3 for T > 200 K (|Burton et al.l 
Il992f) . The S(0) and S(2) lines are from para-H2, while the S(l) and S(3) lines are from ortho-H2. Hence, the line ratios 
of the same "species" are independent of the ortho-to-para ratio while the ratio of two lines adjacent in wavelength is 
sensitive to the ortho-to-para ratio and to the temperature. Unfortunately we have to use a mixture of ortho and para 
lines for observational reasons: a reasonable comparison can only be done for lines observed with the same IRS module 
(slit width). Hence, our analysis of the 'lores' maps will focus on the S(2) and S(3) lines, and the 'hires' spectra on 
the S(l) and S(2) lines. 

In our calculations we assume that the observed line emission is optically thin, the critical densities of the 
observed lines are < 10'^ cm^'^ (this is certainly given for the S(0) through S(2) lines), and that the populated levels 
are in LTE with an ortho-to-para ratio of 3. The S-notation refers to transitions with A J — Jupper — Jiower = +2, and 
with gj = gs{2J + 1) we derive the quantum numbers listed in Table fTOl gs is the nuclear statistical weight factor, 
which is 1 for para-H2 and 3 for ortho-H2. 

In order to derive the excitation temperature T^x we apply the following calculation: the column density iV^ for a given 
transition i is related to the total column density Ntot via the Boltzmann distribution Ni = gjNtot •exp[— E'^/fcTj/Qij^, 
where Qh2 is the temperature-dependent partition function of H2 , with 

Qif, =^{2J + l)g,e-W = 
J 

509.0 1681.7 3474.6 \ / 169.8 1015.1 2503.8 4586.7 

' T + 5e r- +9e ~ + 13e ~ + . . . ) + 3 ■ ( 3e r- + 7e ~ + lie ~ + 15e ~ + . 



(Al) 

For a temperature range between 250 — 400 K the partition function typically varies from 6 to 10. Hence, the 
ratio between column densities for a given transition is 

— = — e >=T . 

N2 32 

The line intensity Ii in [Wm^^] is related to the column density via Ni ~ ^f^^, and the Einstein coefficient Ai 
describes the spontaneous emission probability. Thus the excitation temperature Tex can be calculated as 

„ _ _ E2- El 

J- ax — 



A: In 



32 Al Ji A2 
gi A2 I2 Al 



We emphasize that the mid-infrared H2 lines are only sensitive to the warm (> 100 K) gas, which is only a small 
fraction of the total H2 gas (see Section 221) • In order to derive the mass of warm H2 we use the following calculation: 
for an ortho-to-para ratio of three the total mass of molecular hydrogen is — ^Mortho — wj/j A''//^, where is 
the molecular mass of hydrogen, and the total number of molecules is 

Nh2 



hcAuixj 



xj describes the fractional population in the upper level as xj = Nj/Ntot = gj^^V [^Ej/ht] /Qh2- Combining the 
above equations we get for the warm H2 mass: 

ATrmH2D^^EQH2 
MH2 = ^wj — • 

hcAuigje 
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